N PS  ARCHIVE 
1962 
FLANAGAN,  D. 


AUTOMATIC  CONTROL  OP  PROPELLER  PITCH 

fcy 

LIEUTENANT  DAVID  B.  FLANAGAN 
U.S.  COAST  GUARD 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE 
REQUIREMENTS  FOR  THE  DEGREE  OF  NAVAL  ENGINEER 

and 
FOR  THE  DEGREE  OF  MASTER  OF  SCIENCE  IN 
NAVAL  ARCHITECTURE  AND  MARINE  ENGINEERING 

at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
20  May  1962 


Thesis  Supervisor:   S.  Curtis  Powell,  LOTT.  ING. 
Title:   Associate  Professor  of  Marine  Engineering 


Library 

U.  S.  Naval  Postgraduate  School 

Mon 


AUTOMATIC  CONTROL  OP  PROPELLER  PITCH 

by 

DAVID  B.  FLANAGAN,  U.S.  GOAST  GUARD 
B.S.,  U.S.  COAST  GUARD  ACADEMY 
(1955) 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE 
REQUIREMENTS  FOR  THE  DEGREES  OF  NAVAL  ENGINEER 

and 
MASTER  OF  SCIENCE 
IN  NAVAL  ARCHITECTURE  AND  MARINE  ENGINEERING 

at  the 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
June,  1962 


ETbrary 

U.  '      '  ""  ' 

AUTOMATIC  CONTROL  OP  PROPELLER  PITCH,  by  DAVID  B.  FLANAGAN 
Submitted  to  the  Department  of  Naval  Architecture  and 
Marine  Engineering  on  19  May  1962  in  partial  fulfillment 
of  the  requirements  for  the  degree  of  Naval  Engineer  and 
the  degree  of  Master  of  Science  in  Naval  Architecture  and 
Marine  Engineering. 


ABSTRACT 


The  object  of  this  study  was  to  determine  the  feasi- 
bility of  automatically  controlling  propeller  pitch  to 
reduce  RPM  fluctuations  caused  by  variations  in  propeller 
speed  of  advance.   The  Investigation  was  made  ignoring  the 
reduced  frequency  effect  and  power  plant  governing. 

A  Victory  Ship  (AP2)  was  chosen  as  a  prototype  and  a 
feedback  system  controlling  the  pitch  of  a  controllable 
pitch  propeller  was  developed.   A  non-linear  analogue 
computer  program  approximating  the  system  and  the  propeller 
characteristics  is  presented.   Electronic  difficulties 
prevented  obtaining  total  system  response  from  this  program. 
The  computer  study  Indicated  that  a  linear  analysis  would 
be  useful. 

The  linear  analysis  indicated  that,  for  the  prototype, 
in  order  to  achieve  a  reduction  in  RPM  fluctuation  of  10$ 
of  that  present  in  the  fixed  pitch  case,  a  hydralic  pump 
consuming  over  300  HP  would  be  needed.   Systems  utilizing 
RPM,  shaft  torque,  and  thrust  are  analyzed.   Thrust  sensing 
provided  the  poorest  response  and  in  fact  increased  the  RPM 
fluctuation  for  lower  frequencies. 

RPM  governing  by  automatic  control  of  propeller  pitch 
was  found  to  be  inadvisable  unless  major  design  changes 
are  made  to  reduce  the  spindle  torque  of  presently  available 
controllable  pitch  propellers. 


Thesis  Supervisor:    S.  Curtis  Powell 

Title:  Associate  Professor  of  Marine  Engineering 
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I. 

INTRODUCTION 

The  propulsion  systems  of  ships  operating  In  a  seaway- 
undergo  speed  fluctuations  caused  by  variations  in  propeller 
load.   The  shaft  torque  produced  by  the  propeller  is  a 
function  of  the  flow  velocity  across  the  propeller  disc 
area.  As  the  ship  rolls,  pitches,  and  encounters  waves 
this  flow  velocity  fluctuates  causing  a  fluctuation  in 
shaft  torque  which  the  propulsion  system  must  overcome. 
The  propulsion  system  responds  to  the  fluctuating  load  by 
allowing  fluctuations  in  RPM.   The  amount  of  RPM  fluctuation 
present  depends  on  the  system  inertia,  damping,  and  the 
effectiveness  of  the  governing,  if  provided,  of  the  power 
plant.   High  speed  turbine  installations  with  large 
reduction  gears  experience  less  RPM  fluctuation  difficulties 
than  do  direct  drive  diesels.* 

The  largest  variation  of  flow  velocity  averaged  over 
the  propeller  disc  area,  henceforth  called  speed  of  advance, 
will  probably  occur  at  resonant  pitching.   Model  tests 


*  Evans,  Harvey  J.   Conversations  before  seminar  group 
of  Department  of  Naval  Architecture  and  Marine  Engineer- 
ing, M.I.T.  24  April  1962. 


indicate  that  at  this  condition  the  variation  of  speed  of 
advance  could  be  as  high  as  0.43  of  the  steady  state 
value.    For  a  geared  turbine  ship  such  as  the  Victory 
Ship  (AP2)  with  a  fixed  pitch  propeller  a  fluctuation  of 
about  6  RPM  will  result  at  the  frequency  of  resonant 
pitching, which  is  about  0.03  cycles  per  second.   Calcu- 
lations leading  to  these  figures  are  given  in  Appendix  E. 
For  systems  with  less  effective  inertia  and  damping  the 
fluctuations  would  be  correspondingly  higher. 

Ships  propellers  are  designed  to  operate  efficiently 
at  a  designed  speed  of  advance.   Their  efficiency  is 
reduced  when  subjected  to  a  constantly  varying  speed  of 
advance.   The  propulsion  system  could  deliver  full  power 
when  the  ship  was  operating  in  a  seaway  if  the  load  could 
be  made  sufficiently  constant. 

A  controllable  pitch  propeller  has  the  ability  to 
alter  the  propeller  pitch  at  any  time.   If  the  pitch  were 
altered  to  counteract  the  effect  of  the  changing  speed  of 
advance,  a  constant  torque  load  would  be  provided  to  the 
power  plant.   The  efficiency  of  the  power  plant  and 
propeller  might  be  increased  and  the  entire  system  could 


#*   Burtis,  Evenson  M.  and  J.D.  O'Connell.   Measurement 
of  Ship  Model  Wake  Variation  in  Rough  Water,  Naval 
Engineers  Thesis,  Department  of  Naval  Architecture  and 
Marine  Engineering,  M.I.T.  Cambridge,  1959. 
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be  operated  in  heavy  seas  at  high  power.   A  literature 
survey  has  revealed  that  no  reported  investigation  of 
this  possibility  has  been  made.   Extensive  use  of  one 
governor  systems  controlling  RPM  by  altering  both  pitch 
and  throttle  setting  has  been  made  by  the  aircraft 
industry  in  conjunction  with  turboprop  aircraft. 

The  purpose  of  this  paper  was  to  investigate  the 
possibility  of  governing  propulsion  system  RPM  by  auto- 
matic control  of  propeller  pitch.   Systems  sensing  RPM, 
torque,  and  propeller  thrust  were  investigated.   System 
speed  governing  by  controlling  the  load  may  or  may  not 
be  more  advantageous  than  conventional  governing  of  the 
power  plant.   The  regulation  of  some  variable  other  than 
speed  might  also  provide  a  more  optimum  overall  operation. 
This  investigation,  however,  examines  the  feasibility  of 
speed  governing  through  propeller  load  control  only. 
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II. 

GENERAL  PROCEDURE 

Several  major  analytical  assumptions  were  made  at 
the  outset.   Unsteady  flow  was  treated  using  the  quasi- 
steady  approximation  and  the  reduced  frequency  effects 
were  neglected.   The  speed  of  advance  was  assumed  to  be 
a  steady  state  value  with  a  sinusoidal  fluctuation  super- 
imposed thereon.   The  control  systems  Investigated 
contained  a  simple  one  loop  feedback  with  no  special 
compensating  networks.   The  power  plant  was  assumed  to 
deliver  constant  torque  and  no  model  of  the  governor  was 
developed. 

The  standard  method  of  presenting  propeller  character- 
istics envolves  a  plot  of  non-dimensional  thrust  and  torque 
coefficients  as  functions  of  advance  ratio,  J,  and  pitch 
ratio   .   These  plots  are  very  nearly  families  of  straight 
lines.   For  the  non-linear  study  the  propeller  curves 
were  assumed  to  be  families  of  straight  lines.   The  slopes 
and  values  of  the  actual  and  the  assumed  curves  are  the 
same  at  the  operating  point.   The  accuracy  of  the  approxi- 
mation decreases  as  conditions  depart  from  the  operating 
point.   The  approximations  lead  to  expressions  (6)  and  (13) 


for  propeller  torque  and  thrust  respectively  as  derived 

in  Appendix  A. 

Qp  =  KiPOO2  -  K2Y060   -  K5  0J2       .     (6) 
T  =  K4p^2  -  K5V0^  -  K6  CaJ2         .    (13) 

The  K's  art  constants  for  a  given  propeller  and 
depend  on  both  the  propeller  curves  and  the  dimensions 
of* the  propeller. 

For  the  linear  analysis,  the  propeller  torque  and 
thrust  were  assumed  to  be  linear  functions  of  speed  of 
advance,  pitch,  and  RPM.  This  is  the  customary  incre- 
mental linear  analysis  and  is  valid  only  for  small 
variations  of  the  variables  about  the  operating  point. 
The  equations  for  incremental  variations  in  propeller 
torque  and  thrust  are  derived  in  Appendix  £  and  are 

Qp  =  Kx*vi  +  K*/f>'  +  K*c/j  (81) 

and 

T  =  K*Vq  +  K*/f>'+  K*o>'       .         (100) 

*, 

The  K's  are  found  by  assuming  a  small  change  in 

the  corresponding  variable  (J  0.1),  and  while  keeping 
all  other  independent  variables  constant,  determining 
from  the  propeller  curves  the  corresponding  change  in 
the  dependent  variable.   The  K*  then  becomes  the 
dependent  variable  change  divided  by  0.2   The  Troost 
Standard  Series  for  Single  Screw  Merchant  Vessels  (B4-55) 
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ft 

propeller  curves  were  used  throughout. 

No  test  data  providing  blade  spindle  torque  directly 
is  available.   The  torque  about  the  blade  spindle  axis 
was  assumed  to  be  composed  of  hydrodynamic,  sliding 
friction,  viscous  friction,  inertia,  and  centrifugal 
force  components.  As  with  the  propeller  torque  and 
thrust  coefficients,  the  hydrodynamic  component  of  blade 
torque  can  be  represented  by  a  family  of  blade  torque 
coefficient  curves,  Kfi,  where  advance  ratio  and  pitch 
ratio  are  the  independent  variables.   These  curves  were 
also  assumed  to  be  families  of  straight  lines  and  the 
resulting  expression  for  hydrodynamic  blade  torque  as 
derived  in  Appendix  A  is 

QBH  =  K7PW2  -  K8V0u)  -  K9o02   .         (21) 

Reference  (11)  contains  data  onrche  control  rod 
forces  for  a  controllable  pitch  propeller  Installation 
on  a  tug /(YTB- 502 ),  while  underway  at  various  speeds. 
This  data  also  includes  details  of  the  pitch  setting 
mechanism,  ship's  speed,  centrifugal  force  components  as 
functions  of  RPM  and  pitch,  and  other  useful  information. 
An  approximate  wake  fraction  of  0.164  was  assumed  based 
on  the  block  coefficient  of  the  tug.   Since  the  tests  were 


*  Lewis,  F.M.  "Propeller  Coefficients  and  the  Powering 
of  Ships,"  Society  of  Naval  Architects  and  Marine 
Engineers  Paper,  Nov.  1951. 
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made  under  steady  conditions,  the  only  forces  assumed 
acting  were  hydrodynamic  and  centrifugal.   After 
subtracting  the  centrifugal  component,  the  control 
rod  force  was  converted  to  a  spindle  torque,  non- 
dlmensionalized,  and  finally  plotted  as  blade  torque 
coefficient  versus  advance  ratio  and  pitch  ratio.   The 
results  are  shown  in  Figure  B-3  below  and  the  calcu- 
lations leading  to  this  plot  are  given  in  Appendix  B. 


as 


/S-i 
A2- 
/.!■ 
hO 

•  9- 
.0- 

.*- 

.5- 

3- 
2- 
»/- 


Torque  temdiag  to  rotate  blade  to 
ahead  pitch 


.  f 


T 


J    nD 
FIGURE  B-3 
Derived  Kfi  Curves   for  YTB-502 
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As  before  stated,  the  Kg  curves  were  assumed  linear. 
The  hydrodynamic  spindle  torque  is  a  strong  function  of 
the  location  of  the  spindle  axis  relative  to  the  hydro- 
dynamic  center  of  pressure.   It  was  assumed  that  the 
effect  of  moving  the  spindle  axis  would  be  to  shift  the 
KB  values  up  or  down  the  family  of  Kg  curves.   Once  the 
operating  point  of  the  prototype  assumed  for  this  study 
was  chosen,  the  operating  point  of  the  KB  curves  given 
in  Figure  B-3  was  assumed  moved  to  the  prototypes 
simply  by  shifting  the  Kg  values  corresponding  to  the 
tugs  operating  point  to  the  approximate  position  of  the 
prototype  operating  point.   This  can  he  Justified  only 
in  that  such  a  propeller  could  be  made.   The  slope  and 
spacing  of  the  Kg  curves  was  assumed  to  remain  the  same. 

The  centrifugal  force  component  of  blade  torque  as 
derived  by  Boswell  and  others  is 


r 


Qbc  =  "  V2n 


tip       A   tip 


A   .      \  ^ *max  -  ^min ) 
*xydr  +   \ L 

Hub      J  root 


This  expression  would  be  difficult  to  compute  on  an 


*  Robert  J.  Boswell,  A  Method  of  Calculating  the  Spindle 
Torque  of  a  Controllable-Pitch  Propeller  at  Design 
Conditions.   David  Taylor  Model  Basis  Report  No.  1529, 
Aug.  1961,  pp.  31-34. 
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analogue  computer  and  a  simplification  given  in  Appendix 
A  leads  to  expression  (35)  where  the  constants  K1Q  and  K-^ 
are  functions  of  propeller  geometry  alone. 


QBC  -  "  Kin00  2  "  K;rL^2P  • 


(35) 


The  sliding  friction  component  of  spindle  torque  was 
obtained  from  altering  the  usual  expression  for  sliding 
friction  in  a  Journal  bearing  to  allow  for  a  moment 
tending  to  rotate  the  spindle  shaft  about  a  line  perpen- 
dicular to  the  shaft  centerline.   The  resulting  expression 


is 


Qbs  =  -    S 


(38) 


The  double  sign  in  expression  (38)  Indicates  that  the 
direction  of  the  torque  must  be  taken  either  (+)  or  (-) 
depending  on  whether  the  direction  of  pitch  change  is 
toward  lower  or  higher  ahead  pitch  respectively. 

The  remaining  components  of  blade  spindle  torque  are 

obtained  by  direct  application  of  their  respective 

definitions.   The  final  expression  for  blade  torque  is 

hydrodynamic 


QB  =  'KypfcO2   -  K8V0  CO  -  K9  6/J  2 


sliding  friction       viscous   friction 


centrifugal 


-     K10U^   -  K1;LcAJ2p        (38) 
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There  are  numerous  variations  of  hub  mechanisms  in 
use.   All  are  required  to  convert  motion  in  the  direction 
of  the  propeller  shaft  axis  into  rotation  of  the  blades 
about  their  spindle  axes.   As  shown  in  Appendix  A,  each 
of  the  four  mechanisms  most  frequently  used  result  in 
approximately  the  same  amount  of  blade  rotation  for  the 
same  input  of  linear  motion  along  the  propeller  axis  up 
to  a  pitch  ratio  of  1.5.   The  same  is  true  for  the  force- 
to-torque  transmission.   For  this  reason  the  simplest 
mathematical  relation  of  the  four  was  used  throughout. 
The  results  are  approximately  applicable  to  all  the 
mechanisms.   The  expressions  used  are  as  follows: 

=  lF  (39) 


and 


FCQ  =hB       •  <41> 


R 

There  are  two  general  types  of  controllable  pitch 
propeller  mechanisms.   The  first,  designated  herein  as 
Design  I,  utilizes  a  control  rod  extending  through  the 
propeller  shaft  and  is  driven  by  a  hydralic  piston  with- 
in the  ship.   The  second,  designated  herein  as  Design  II, 
is  driven  by  a  power  piston  housed  within  the  propeller 
hub.   The  power  piston  of  Design  II  receives  oil  from 
a  pump  within  the  ship  through  oil  passages  in  the 
propeller  shaft.   Design  II,  as  presently  available, 
includes  a  safety  spring  within  the  hub  that  makes  the 
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propeller  "rail-safe"  by  forcing  the  blades  into  ahead 
full  pitch  in  the  event  of  a  failure  in  hydralic  oil 
pressure. 

The  equation  of  motion  of  both  these  designs  is 
obtained  by  equating  the  oil  flow  from  the  pump  with 
the  sum  of  the  flows  into  the"  mechanism  which  include 
flow  from  piston  motion,  flow  from  leakage,  and  flow 
due  to  compression  of  the  fluid  and  expansion  of  piping. 
The  final  equation  of  both  systems  are  of  the  same  form 
with  the  exception  of  additional  terms  in  the  Design  II 
expression  that  are  generated  by  a  safety  spring  force. 
Design  I  was  used  for  all  further  study.   Schematics  of 
both  designs  are  given  as  Figures  A-7  and  A-8  in  Appendix 
A.  .  The  equation  of  motion  of  Design  I  as  derived  in 
Appendix  A  is 

Input     Ram  Flow       Leakage  flow 

K7p  tp  2  -  KgVJO 

R        R 


Kpx  =  p  i2^  +  *k 

D     AR 


L 

-  l9±±l     J   (fR  ♦  **£-*)   ♦-*  t£*   W> 

R  D*      TTD     R 


+  L.  (£iM  +  McR)  -  Kioa2  ~  Kll  (£  2P 
-n-D    e^         R        R 


Compliance  Flow 
+  ^2.  )   -1  [  /  K?p  OJ  2   -  KsVqUJ  -  Koco2 
AH   «   V    e        ~T~       R 
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bR   /     R         R 

]L  (<££  V'R*   )+  V   (JiJ  +  McR^l    (42) 
D    R  1TD   -J—     J 

The  Victory  Ship  (AP2)  was  chosen  as  a  prototype  to 
determine  values  of  the  various  parameters.   This  is  a  6000  HP 
geared  turbine  merchant  ship  for  which  considerable  data  was 
available.   The  detailed  determination  of  the  various 
parameters  present  in  equation  (42)  above  is  carried  out  in 
Appendix  B.   The  propeller  blade  was  assumed  to  be  an 
elliptical  flat  plate  for  the  determination  of  blade  inertia 
and  centrifugal  force  components  of  blade  spindle  torque. 
The  dimensions  of  the  hub  were  obtained  by  scaling  up  the 
design  given  by  Rupp  in  Reference  11. 

The  three  remaining  parameters  to  fix  were  area  of  the 
ram,  the  size  of  the  pump,  and  the  hydrallc  pressure.   Since 
the  initial  cost  and  maintenance  on  systems  using  more  than 
1000  psi  would  be  high,  this  value  was  chosen  as  the 
maximum  hydralic  pressure.   The  maximum  diameter  of  power 
piston  that  could  be  accommodated  within  the  propeller  hub 
has  a  radius  of  15.75  inches. 

Although  Design  I  has  the  power  piston  within  the  ship, 
this  was  assumed  the  largest  piston  size  that  would  be  fitted. 
Calculations  given  in  Appendix  B  show  that  with  the  assumed 
hydralic  pressure  of  1000  psi  and  this  ram  area  there  is  still 
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a  likelihood  of  stalling  the  mechanism  (by  opening  relief 
valves)  even  with  this  large  siz.e  power  piston.   The 

maximum  combinations  of  blade  centrifugal,  hydrodynamic , 

p 
and  sliding  friction  when  acting  together  require  1000  in 

of  piston  area.   The  assumed  piston  area  is  775  in2. 

Larger  piston  sizes  would  require  larger  pump  capacities 

which  are  excessive  for  this  piston  area. 

As  shown  in  Appendix  E,  the  time  constant  of  the  pitch 

control  mechanism  is 

^P  =  — B .  (60) 

KpTTpM 

The  effectiveness  of  the  system  in  regulating  RPM 
fluctuations  requires  that  this  time  constant  be  small. 
The  larger  pump  capacities  correspond  to  larger  Kp.   Assuming 
that  the  marine  industry  would  not  accept  a  device  consuming 
more  than  5%   of  the  installed  main  plant  horsepower,  a 
maximum  pump  size  of  300  HP  was  assumed.   Kp  was  then 
determined  from  manufacturers  data. 

The  viscous  friction  forces  were  obtained  from  an 
analysis  of  transient  spin  test  data  taken  by  Rupp.   This 
data  consists  of  a  plot  of  control  rod  forces  versus  pitch 
for  the  operating  conditions  of  going  from  full  ahead  to 
full  astern  and  vice-versa.   This  data  is  reproduced  as 
Figure  B-7  below.   The  thrust  forces  causing  high  values 
of  sliding  friction  were  assumed  zero  for  this  test.   The 
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FIGURE  B-7 

YTB-502  Control  Rod  Force  for  Spin  Test 
In  Air  From  Data  by  Rupp 

long  region  of  both  transient  curves  that  are  parallel 

to  the  steady  state  condition  are  assumed  to  correspond 

to  the  condition  of  maximum  pump  output  and  constant  piston 

speed.   The  deviation  from  the  steady  state  curve  in  both 

cases  is  approximately  3  x  10  lbs.   This  control  rod 

force  was  converted  to  spindle  torque  and  scaled  up  by  the 

ratio  of  propeller  diameter  squared  since  viscous  friction 

is  proportional  to  surface  area.   The  speed  of  rotation 

was  assumed  constant  and  equal  to  the  total  angle  divided 

by  total  time. 


j 
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The  substitution  of  the  various  parameters  into 
equation  (42)  revealed  that  several  of  the  terms  would 
be  very  small  in  comparison  to  others.   The  coefficients 
of  each  independent  variable  was  examined  and  those 
terms  having  coefficients  of  less  than  1%   of  those  re- 
maining for  the  same  variable  were  eliminated.   A 
further  simplification  was  made  by  setting  to  zero  all 
terms  including  the  first  derivative  of  shaft  speed. 
This  assumption  is  based  on  the  knowledge  that  propeller 
RPM  does  not  vary  to  a  large  extent  even  without  regu- 
lation and  should  vary  less  with  regulation.   The 
resulting  simplification  of  equation  (46)  is 
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III. 

PROCEDURE  AND  RESULTS  OP 
THE  NON-LINEAR  ANALOGUE  COMPUTER  STUDY 

A.  Procedure 

Equations  (6),  (13) »  and  (65)  given  above  envolve 
non-linearities  in  the  form  of  products  of  independent 
variables  and  a  discontinuous  sliding  friction  term.   If 
saturation  of  the  pump  is  also  allowed,  another  discon- 
tinuity is  introduced  which  would  limit  the  value  of  pump 
output,  kpx,  to  some  finite  value.  All  of  these  non- 
linearities  can  be  programed  on  the  analogue  computer 

that  was  available.  The  program  includes  four  multipliers 

2     2  ^oJJ      * 
to  generate^  »  p  60  ,Aand  v060  .   It  also  includes  a 

switching  mechanism  controlled  by  the  derivative  of  pitch 

to  approximate  sliding  friction. 

The  program  was  arranged  so  that  shaft  torque  and 
thrust  were  generated  in  order  to  test  the  final  system 
sensing  these  variables. 

The  program  is  shown  in  Figure  C-2  below.   The  para- 
meters determined  in  Appendix  B  were  then  inserted  into 
the  program  and  voltage  scales  were  assumed  for  each 
variable.   The  gains  of  various  loops  were  then  adjusted 
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to  provide  signals  that  would  not  saturate  the  various 
amplifiers  and  multipliers  of  the  computer.   The  details 
of  this  scaling  is  given  in  Appendix  C. 

The  input  to  the  system  was  taken  to  be  the  voltage 
analogue  corresponding  to  speed  of  advance,  VQ,  and 
consisted  of  a  steady  positive  value  with  a  sinusoidal 
voltage  of  amplitude  &V0  provided  by  a  signal  generator. 
The  output  consisted  of  a  voltage  analogue  for  shaft 
speed,  to  ,  and  both  input  and  output  were  recorded  on  a 
two  channel  Sanborn  recorder.   It  was  hoped  that  data  over 
a  range  of  frequencies  could  be  obtained  in  this  manner 
for  the  fixed  pitch  case  as  well  as  the  system  operating 
in  closed  loop  sensing  shaft  speed,  torque,  and  thrust. 

B.  Results  of  the  Non-Linear  Study 

The  frequency  response  of  the  fixed  pitch  system  was 
obtained.   The  data  is  recorded  in  Appendix  C  and  is 
shown  graphically  as  a  plot  of  the  magnitude  of  the  ratio 
of  output  to  input  in  decilogs  in  Figure  I  below.   The 
strong  stabilizing  effect  of  sliding  friction  was  very 
apparent.   When  the  signal  corresponding  to  sliding 
friction  feedback  was  removed  with  a  sinusoidal  Av0 
input  of  amplitude  6  ft/sec  the  second  derivative  of 
pitch,  p%  increased  by  a  factor  of  80.   The  feed  forward 
path  containing  Vo^aj   had  no  effect  on  p   even  for  very 
large  values  of   kVQ. 
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Values  generated  by  the  computer  for  propeller 
torque,  thrust,  and  blade  torque  were  within  5$  of  the 
calculated  values.   The  time  constant  of  the  system  with 
fixed  pitch  was  observed  by  injecting  a  step  in  engine 
torque  and  recording  the  decay  in  shaft  speed  variation. 
The  time  constant  was  observed  to  be  approximately  3.28 
seconds.   The  decay  to  steady  state  required  about  10 
seconds  (approximately  3  time  constants). 

Electronic  difficulties,  the  details  of  which  are 
given  in  Appendix  D,  prohibited  obtaining  data  for  the 
control  system  operating  in  closed  loop. 

In  all  cases  the  node  containing  p  was  very  small 
and  consisted  primarily  of  computer  generated  noise. 
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IV. 
PROCEDURE  AND  RESULTS  OF  THE  LINEAR  STUDY 

A.  Procedure 

The  results  of  the  non-linear  study  Indicated  that 
valid  further  simplifications  could  be  made  in  equation 
(46).   Eliminating  the  terms  containing  Vtu    and  p  as 
suggested  by  Chapter  III  above,  equation  (46)  reduces  to 


RAr  \   Kpx  -  5i 


(k7-ku)  E±L       (67) 

R 


-  K8Vo^   -  (K9+K10)  ^2    (f^hllAj.   )~R 
R  "R  bR     J) 

Equation  (67)  includes  ram  flow  and  leakage  terms 
only.   Since  the  linear  analysis  is  valid  only  for  small 
variations  around  the  operating  point,  the  leakage  was 
assumed  constant  and  is  equivalent  to  an  effective 
reduction  in  pump  output  Kp#.   As  shown  in  Appendix  E, 
equation  (67)  eventually  reduces  to  a  formula  whose 
frequency  domain  transfer  function  is 

y»(s)   V+1 
where  y'(s)  =  Laplace  transform  of  the  incremental 


mechanical  input  to  the  pitch  control  mechanism. 

Using  the  linear  approximations  for  propeller 
response  before  mentioned,  and  assuming  the  power  plant- 
propeller  system  to  be  a  first  order  system  containing 
Inertia  and  damping,  the  complete  systems  for  sensing 
RPM,  torque,  and  thrust  were  drawn  in  block  diagram 
form  as  shown  in  Figure  II. 

The  gains  G^  and  Gq  for  the  shaft  speed  and  shaft 
torque  sensing  systems  respectively  were  chosen  so  that 
a  damping  ratio  of  0.4  was  assumed  by  the  closed  loop 
system.  Higher  values  of  loop  gain  reduce  the  low 
frequency  response  but  cause  higher  peaking  or  ampli- 
fication at  a  higher  frequency.   This  peaking  is  apparent 
at  about  2  radians  per  second  as  shown  in  Figure  E-6.  The 
gain  for  thrust,  GT,  was  taken  to  be  a  very  low  value 
only  to  demonstrate  that  any  positive  value  of  Gm  is  not 
beneficial  in  regulating  RPM  a.%   the  lower  frequencies. 

B.   Results 

System  gain  is  defined  as  the  ratio  of  shaft  speed 
variation  for  a  given  variation  in  speed  of  advance,  W  . 
The  system  gain  for  each  of  the  three  linear  systems 
pictured  in  Figure  II  was  obtained  and  plotted  by  the 
standard  methods. 

The  results  are  shown  in  Figure  E-6  where  system  gain 
is  plotted  in  decilogs. 
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V. 
DISCUSSIOM  OF  RESULTS 

As  shown  in  Figure  E-6,  the  fixed  pitch  response 
determined  on  the  analogue  computer  was  in  good  agreement 
with  the  linear  analysis  for  frequencies  below  0.10  rad/sec. 
It  is  of  interest  that  the  input  levels  used  to  obtain  this 
response  curve  are  extremely  high.   The  minimum  value  for 
the  fluctuating  component  of  the  speed  of  advance  is  three 
volts  which  corresponds  to  thirty  feet  per  second.   The 
steady  state  value  of  speed  of  advance  is  only  16.7  ft/sec. 
The  maximum  value  of  input  reached  40  volts  corresponding 
to  400  ft/sec.   These  high  values  of  input  were  necessary 
to  obtain  an  output  above  the  noise  level.   The  results 
therefore  indicate  that  the  linear  analysis  should  be 
reasonably  valid  even  for  moderate  fluctuations  around  the 
operating  point. 

Had  it  been  known  that  the  terms  containing  p  and  V  in 
equation  (65)  would  be  insignificant,  an  analogue  eliminating 
these  paths  would  have  been  developed.   As  explained  in 
Appendix  D,  hum  in  the  *p  path  made  the  obtaining  of  data 
for  the  closed  loop  system  impossible.   The  simpler  analogue 
would  not  have  this  low  signal  level  path  and  might  have  been 


more  successful.   Unfortunately  time  was  not  available 
to  reprogram. 

The  simplifications  leading  to  equation  (65)  as 
shown  in  Appendix  C,  call  for  ignoring  the  term 

SL  (fR  *  fr?±*   g  ) 
dt  bR 

This  friction  term  represents  an  oil  flow  and  since 
it  is  discontinuous  at  p  =  o,  corresponds  to  an  impulse  or 
hydralic  shock.  This  hydralic  shock  would  very  likely  be 
damaging  to  relief  valves  and  other  parts  of  the  system. 
The  introduction  of  an  accumulator  to  cushion  this  shock 
would  introduce  additional  compliance  into  the  system,  but 
as  shown  in  Appendices  C  and  D  compliance  is  relatively 
unimportant.   There  are  no  compliance  terms  remaining  in 
the  final  linear  analysis. 

The  simplifications  of  expression  (46)  given  in 
Appendix  C  show  that  inertia  loading  is  relatively  unimpor- 
tant. Design  I  has  a  larger  Inertia  loading  because  of  the 
relatively  larger  mass  of  the  control  rod  compared  to  the 
smaller  power  piston  of  Design  II.   The  area  of  the  high 
pressure  port  delivering  full  flow  to  the  power  piston  of 
Design  II  would  have  to  be  9.1  square  inches  to  match  the 
output  of  the  variable  displacement  pump  at  the  same  pressure. 
Since  inertia  loading  is  relatively  unimportant  and  since 
providing  a  large  oil  port  area  would  be  a  difficult  engineer- 
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ing  problem,  there  is  little  advantage  to  Design  II  for 
this  application. 

The  transfer  functions  used  in  comparing  the  various 
systems  are  stable  for  all  frequencies.  These  transfer 
functions  are  based  on  numerous  assumptions  that  neglect 
many  of  the  lags  and  storage  elements  actually  present. 
They  also  do  not  consider  nonlinearities  such  as  sliding 
friction,  st^iction,  saturation,  and  the  nonlinearities 
present  in  the  propeller  characteristics  themselves. 
There  is  therefore  a  possibility  that  the  actual  hardware 
might  exhibit  natural  frequencies  corresponding  to 
resonances  in  the  remainder  of  the  propulsion  system. 
The  lowest  frequency  of  torsional  vibration  resonance  for 
the  Victory  Ship  (AP2),  as  measured  and  reported  in 
Reference  16,  is  28  cycles  per  second  or  175  radians  per 
second. 

With  the  restrictions  of  a  maximum  pressure  of  3,000 
psi.  hydralic  pressure,  a  maximum  pump  power  of  300  HP, 
and  the  approximate  dimensions  of  present  day  controllable 
pitch  propellers,  the  improvement  in  shaft  RPM  regulation 
over  a  fixed  pitch  propeller  installation  using  automatic 
pitch  control  corresponds  to  a  reduction  in  RPM  fluctuations 
by  a  factor  of  nine.   The  improvement  in  propeller  and 
power  plant  efficiency  would  not  be  great  enough  for  this 
small  an  increase  in  RPM  regulation  to  warrant  the 
installation  of  such  a  large  and  complex  system. 
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Several  improvements  in  the  system  analyzed  might 
reduce  the  size  of  the  equipment  necessary  or  improve 
the  system  response.   The  raising  of  the  hydralic 
pressure  would  reduce  the  size  of  the  piston  necessary 
and  thereby  reduce  the  capacity  of  flow  necessary  from 
the  pump.   This  would  improve  the  system  by  reducing 
the  important  pump-piston  time  lag.   Increasing  the 
horsepower  of  the  pump  also  improves  the  system  by 
decreasing  the  pump-piston  time  lag.   Compensation  in  the 
form  of  lead  or  lag  networks  in  the  feedback  loop  would 
delay  the  peaking  present  in  the  proposed  systems  and 
allow  a  higher  closed  loop  gain  and  hence  a  higher  degree 
of  regulation. 

It  is  unlikely  that  the  above  mentioned  changes 
would  result  in  improvementsJLn  system  performance  to  any 
large  degree.   If  spindle  torque  could  be  greatly  reduced, 
the  size  of  the  power  piston  could  be  reduced  and  the 
entire  system  could  be  designed  wiTth  better  response 
and  be  made  smaller.   As  shown  by  equation  (38),  spindle 
torque  consists  of  many  components,  each  of  them  large. 
By  increasing  hub  size  the  spindle  bearing  could  be  made 
larger  and  the  lever  arm  of  the  hub  mechanism  increased. 
The  application  of  non-friction  bearings  might  reduce 
this  factor  although  roller  and  ball  bearings  are  presently 
considered  unacceptable  for  naval  service  due  to  shock 
susceptability.   If  the  propeller  were  designed  to  operate 
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at  one  RPM,  the  hydrodynamlc  and  centrifugal  components 
might  be  reduced  by  careful  design.   These  significant 
changes  might  make  automatic  control  of  pitch  to  remove 
fluctuation  in  RPM  feasible. 

As  demonstrated  by  Figure  E-6  there  is  no  advantage 
in  sensing  either  torque  or  thrust  in  lieu  of  RPM  when 
the  variable  being  regulated  is  RPM*  The  Improvement 
in  regulation  offered  by  RPM  sensing  is  9.7  decilogs 
which  corresponds  to  a  reduction  in  fluctuation  by  a 
factor  of  9. 
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VI. 
CONCLUSIONS 

For  geared  turbine  or  geared  diesel  propulsion  plants 
used  by  the  U.S.  maritime  industry,  automatic  control  of 
propeller  pitch  to  reduce  wave  excited  RPM  fluctuations 
would  be  excessively  expensive  in  horsepower,  space,  and 
complexity  for  a  small  gain  in  RPM  regulation. 

Significant  reductions  in  dynamic  blade  spindle  torque 
through  design  changes  would  be  necessary  to  make  this 
scheme  practical. 


VII. 
RECOMMENDATIONS 

Automatic  pitch  control  investigation  requires  a 
knowledge  of  spindle  torque.  Very  little  is  known  in 
this  area.   It  is  recommended  that  future  effort  be  con- 
centrated on  the  determination  of  dynamic  spindle  torque, 

A  more  useful  approach  to  the  problem  of  automatic 
pitch  control  would  include  the  presence  of  engine 
governing  with  a  view  toward  developing  a  system  that, 
with  one  governor  would  receive  a  bridge  signal  and 
control  pitch  and  engine  throttle  setting  to  achieve  the 
thrust  or  ship's  speed  desired  by  the  operator. 


VIII. 
APPENDIX 


APPENDIX  A 
Derivation  of  System  Equations 

1.   The  Power  Plant-Propeller  System: 

Figure  A-1  represents  a  schematic  of  the  propeller- 
power  plant  system. 

The  equation  of  motion  of  this  system  is 

QE  "  QP  ~  £2^=   (JiP  +  J1^  &  J  <1) 

where 

QE   =  developed  engine  torque,  lb-ft 

Qp   =  absorbed  propeller  torque,  lb-ft 

» 

Jip  =  effective  polar  mass  moment  of  inertia  of  the 

propeller  (includes  entrained  water),  lb-ft-sec2 

1 

Jig  =  effective  polar  mass  moment  of  inertia  of  the 

engine  and  gearing  referred  to  the  propeller  shaft, 
lb-ft-sec2 


°£    -   mechanical  damping  in  the  power  plant  and  gearing, 
*  ft-lb-sec. 


Qs   =  shaft  torque,  ft-lbs. 

CO  =  shaft  rotational  speed,  radians  /  second 
cu  =  shaft  rotational  acceleration,  rad/sec  . 
The  torque  delivered  by  the  engine  is 

QE  -£e^  • 
The  shaft  torque  expression  Is 

QS  =  QE  -«£EW"  J1ijk      '  {2) 

Curves  of  Kq  for  most  propellers  are  nearly  linear  over  a 
wide  range  of  adrance  coefficient  and  pitch  ratio.   Figure  A-2 
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FIGURE  A-l 
The  Propeller-Power  Plant  System 
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FIGURE  2 
Assumed  Linear  Torque  Coefficient  (Kq)  Curves 
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represents  an  approximate  general  set  of  Kq  curves  assuming 
that  they  are  completely  linear. 

The  equations  of  these  lines  are 

5=  "V  +  Yq  ♦  Hq 


f  =     YD- J    +    hK~     A    V,.  (3) 


where 

KQ 


=  a£ 

Pn2  JP 


p   =  propeller  pitch,  ft. 

D   =  propeller  diameter,  ft. 

no  q  =  slope  of  the  Kq  curve  at  the  operating  point  of  the 
true  Kq  curve  and  assumed  constant  for  this  analysis, 


non-dimensional 
J   =  propeller  advance  coefficient  =  ^°- 

Vo  =  speed  of  advance  for  the  propeller,  ft/sec 

n   =  propeller  speed,  rev. /sec 

Af 
bn  =  ratio     D    ,  non-dimensional,  assumed  constant 

hQ  =  t/D  intercept  of  the  KQ=0  line. 

P   =  water  mass  density  1 .024 

When  m  q ,  bo,  and  hQ  correspond  to  those  at  an  operating 
point  taken  from  actual  propeller  curves,  the  response  of  the 
propeller  to  changes  around  the  operating  point  will  be 
accurate  within  the  limitations  of  a  steady  flow  analysis.   As 
the  propeller  operates  farther  from  this  point  the  assumption 
of  linearity  becomes  less  valid,  but  is  considered  sufficiently 


1 .   This  unit  for  mass  density  for  water  was  taken  to  be 
consistant  with  Troost's  propeller  curves  used  to 
determine  coefficients  of  equation  (3). 
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accurate  for  this  study. 

Solving  equation  (3)  for  Qp, 

Qp  =  f  92D4  -  ^qVo  P"D4  -  hQ 9  n2D5        #        (4) 
bQ         bQ         bQ 

Since  n  =  —   ,  (5) 

2TT 


equation  (4)  can  be  written  in  the  form 
where 


«P 

=  K-j/pco2   - 

K2Voa)- 

K3 

«j2 

K1 

4TT*bQ 

» 

K2 

=  rnQ  f  D4 
2TTbQ 

» 

(6) 


(7) 


(8) 

and  K,  =  hp  ?lP  (9) 

Substitution  of  equation  (6)  into  equation  (1 )  completes 
the  derivation  of  the  propeller-power  plant  system  equation. 
QE  -  K^foo2   +  K2Vo£l>  +  K36O2  -cC$cO   =  (jip  +  JiE)  u>  (10) 

The  thrust  coefficient  for  most  propellers  can  also  be 
considered  linear  with  respect  to  f/D   and  J.   Again  this  is 
accurate  for  small  changes  around  an  operating  point  and  less 
accurate  for  large  departures  from  this  point.   Figure  A-3 
represents  the  linear  approximation  to  the  actual  KT  curves. 

The  procedure  for  obtaining  thrust  as  a  function  of  Vo,u;, 
and  p  follows  in  the  same  manner  as  used  in  obtaining  the 
similar  expression  for  Qp. 

Prom  the  linearized  KT  curves 

<f    =       KttJ  +  bTKT+h_  „  (11) 

D  Till 
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FIGURE  A-3 
Linear  Approximation  to  KT  Curves 
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FIGURE  A-4 
Linearized  Hydrodynamic  Blade  Torque  Coefficient  Curves 
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Using  the  definition  of  KT  , 

Kr=    jj&r-        ,  (12) 

and  solving  for   I   . 

T  =  KA/paj2  -  K5Vo^)  -  K6  uO2  (13) 

where 

K5  =   ^Pp3     i  (15) 

and  K6  =   ^T^4"     •  (16) 

With  b  =  hSL  »  (17) 

*T)T  =  slope  of  KT  curve, 
and  hT  =  X  intercept  of  Kq  =  0  line. 


2.   Spindle  Torque 

The  equations  of  motion  of  any  servo  system  are  a  function 
of  system  loading.   Loading  in  this  case  consists  of  blade 
spindle  torque.   Torque  about  the  spindle  axis  can  be  viewed 
as  composed  of  hydrodynamic ,  sliding  friction,  viscous  friction, 
inertia,  and  centrifugal  force  components. 

These  components  are  derived  separately.   The  equation 
for  spindle  torque  is 

Qb  =  Qbh  +  Qbs  +  Qbv  +  QBj  +  QBC  (18) 

where 

Qg  =  total  spindle  torque  that  must  be  overcome  by  the 
pitch  changing  mechanism,  lb-ft 

QBH  =  hydrodynamic  spindle  torque,  lb-ft 

QBS  =  sliding  friction  torque,  lb-ft 
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QBv  =  viscous  friction  torque,  lb-ft 

QBJ  =  J1b  ^B*  lb-ft 

QB0  =  centrifugal  force  torque  from  the  mass  of  the 
propeller  metal,  lb-ft. 

Very  few  curves  of  spindle  torque  are  available.   Hydro- 
dynamic  spindle  torque  is  a  strong  function  of  the  spindle 
axis  position  relative  to  the  center  of  hydrodymic 
pressure.   Tests  made  at  M.I.T.  show  that  curves  of  Kg,  where 
KB  is  defined  as 

KB  =  QBH  (19) 

are  very  nearly  linear  functions  of  advance  coefficient  and 

(2) 
pitch  ratio  just  as  is  KA^   As  with  Q  linear  curves  for  KB 

are  assumed  as  shown  in  Figure  A-4,  and  the  expression  for 

QgH  is  derived  based  on  this  assumption. 

For  these  curves  a  spindle  torque  that  tends  to  increase 
pitch  is  considered  a  positive  torque. 

Proceeding  as  before, 

f=  niBJ  +  bB  KBH  +    hB  (20) 


and 


2     _    tr_\m        _    v     /••>/- 


QBH  =  K7fo)^   -  KgVO     -  K^o)*  (21) 

where 

K7  =      P  D4  (22) 

7    *TO 


2 J     Conversation  with  S.    Curtis  Powell,   March  2,   1962, 
M.I.T. 
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K8  =      "iflpD4  ■  (23) 

2TTbB 

and     K9  =       h*    P    D3  .  (24) 

The  spindle  Is  not  a  journal  loaded  In  the  normal 
manner.   Propeller  thrust  acts  at  a  considerable  moment 
arm  tending  to  rotate  the  spinde  about  an  axis  perpen- 
dicular to  its  centerline  as  shown  in  Figure  A-5. 

If  the  bearing  were  normally  loaded  the  friction  torque 
would  be 

fPr  ,   lb-ft        j 
where 

f  =  coefficient  of  journal  friction 

P  =  bearing  load,  lb 

r  =  spindle  radius,  ft. 

Tl  =  total  thrust  on  blade  perpendicular  to  the  spindle 
axis,  lb. 

An  approximate  expression  for  spindle  sliding  friction  is 

QBS  =  fr  ii_M  (  j 


where 


JL   =  lever  arm  from  center  of  thrust  on  the  blade  to  the 
center  of  the  spindle  bearing,  ft. 

b  =  spindle  bearing  length,  ft. 

B  =  number  of  blades 

Sliding  and  static  friction  are  assumed  to  be  the  same. 

Blade  viscous  torque  can  be  expressed  as 

Qbv  =<£  B^BN  (26) 

-4-0- 


HUB 


<w)s 


s^M  i 


Sp»nde   bearing 


FIGURE  A-5 
Schematic   Indicating  Blade  Sliding  Friction   (QBg) 


2ITR,»  TTD 


FIGURE  A- 6 
The  Relationship  Between  Pitch  and  Blade  Angle 
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where 

°£B  =  coefficient  of  viscous  friction,  lb-ft-sec. 
CO    =  spindle  rotational  velocity,  radians/sec. 
Blade  inertia  torque  is 

QBJ  =  JiB  ^B  N       )  (27) 

where 

Jig  =  polar  mass  moment  of  inertia  of  the  blade  about 
the  spindle  axis,  lb-ft-sec2 

•  p 

cjg   =  spindle  rotational  acceleration,  rad/sec   . 

The  expression  for  blade  spindle  torque  as  derived  by- 
Bo  swell  and  others  is  p  x 

*-  hub  root 


where 


O    =  mass  density  of  the  blade  material, 
\p   lb-sec2/ft4 

60  =  propeller  angular  velocity  in  rad/sec, 

A     =  blade  section  area,  ft2 

x3y  =  coordinates  of  the  blade  section  centroid 
measured  in  the  plane  of  the  section  from 
the  blade  spindle  axis,  ft. 

r   =  radius  of  the  section,  ft. 

•^MIN  =  ^e  minimum. moment  of   inertia  of  the  blade 
section,  ft4 


3.   Robert  J.  Boswell,  A  Method  of  Calculating  the  Spindle 
Torque  of  a  Controllable-Pitch  Propeller  at  Design 
Conditions.   David  Taylor  Model  Basin  Report  No.  152y, 
Aug.  1961,  pp.  31-34. 
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•^MAX  =  tile  maximumAmomen-t  of  Inertia  of  the  blade 
section,  ft^ 

oC  =  pitch  angle  of  the  section  from  its  major 
axis,  and  is  a  function  of  both  r  and  p  „ 

Q_0   =  centrifugal  component  of  blade  torque,  positive 
tending  to  rotate  pitch  ahead.,  lb-ft 

Equation  (28)  is  not  readily  adapted  to  analogue 
computation.  Several  simplifying  assumptions  are  necessary 
to  make  it  so. 

Let  0 -q  be   the  pitch  angle  in  radians  at  the  seven 
tenths  radius.  This  angle  corresponds  to  the  "pitch"  of  the 
blade  as  shown  in  Figure  A-6. 

cC  in  equation  (21)  will  be  greater  than  O      near  the 
root  and  less  than  &  B  at  the  tip. 

Assume  cC  =  0  B.  This  assumes  a  constant  pitch  angle 
over  the  blade  length.  Equation  (28)  can  then  be  written: 

-  Qbo  = 

hub  °  root 


^Nf\A/*u)dr  +  p60z,W£B  \tWl^^   (29) 


From  Figure  A-6  is  can  be  seen  that 

TTD  B 

and      .\       6B=  tan-1    ^  j  (3D 

whose  expansion  is 

9  =.±  -I/^3   +  ...    ,  (32) 


f  or  £  < 

TTD 


The  expansion  for  Sin  2  #B  is 
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Sin  26^  =   2(9  ^  -  (g^B)3  +  ...  (33) 

B       B       3  / 

e 

Equating  the  first  terms  of  these  expansions, 

2_f  =  Sin  2  9^  (34) 

TD  B 

This  assumption  will  introduce  an  error  for  this 

component  of  spindle  torque  that  is  sometimes  conservative 

and  sometimes  not  conservative  since  the  total  spindle  torque 

is  a  sum  of  components  with  changing  signs  leading  to  a  sum 

that  also  changes  sign.   The  maximum  error  will  occur  for 

f  /  4°  \ 

large  values  of  —  .     For  the  case  studies,   (  i  1     =1.1 

B  VD  i  MAX 

For  this  case 

2f  =  0.7  05, 
TTD 

and  Sin  2  0  ^  =  0.644  . 
B 

With  these  assumptions  equation  (29)  reduces  to 

Qbc  =  -Kio^2  "  Kif°2  (35) 


where 


K10=         ^p\     A^dr  (36) 

J  hub 


and 


Kll  = 


0tip 
M^DPdX  (I         -  I        )   dr  (37) 


MAX  MIN' 

'hub 


Substituting  equations    (21),    (25),    (26),    (27),    and    (36) 
into   equation    (18)     . 
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Q     =  K7/pO)2  -  KqVoO£>  -  K9602  i  f  ££Jn    +  *CB  ^bn)    + 

b 

JiB  <£B    -  K10602  -  Kn  OJ2^     .  (38) 

Where  Nl  =  number  of  blades 

The  double  sign  before  the  fourth  term,  (SLLka  ),  is  to 

b 

indicate  that  this  term  changes  its  sign  such  that  friction 
torque  resists  spindle  motion.   This  term  will  take  on  the 
same  sign  as  the  viscous  friction  term,  cC^  oj^  . 

-D     J3 

3.   The  Pitch  Setting  Mechanism 

Two  designs  for  the  pitch  setting  mechanism  are  investi- 
gated.  Both  are  hydralic.  Design  I  is  a  pump  controlled 
system  using  a  variable  displacement  pump  and  a  control  rod 
running  through  the  propeller  shaft  and  is  shown  schematically 
in  Figure  A-7.  Design  II  is  a  valve  controlled  system  with  the 
power  piston  within  the  hub.   Design  II  includes  a  safety  spring 
to  make  the  propeller  fail  safe  in  case  of  loss  of  hydralic 
pressure  and  is  shown  in  Figure  A-8. 

The  relationship  between  pitch  and  control  rod  or  piston 
position  depends  on  the  particular  hub  mechanism  involved. 
Several  designs  are  in  use.   Figures  A-8  and  A-7  show  two 
possible  methods.   Figure  A-7  shows  two  additional  methods 
also  in  use. 

For  a  given  hub  diameter  the  parameter  R  in  each  mechanism 
will  be  approximately  the  same.   The  ratio  §2.  for  a  given  pitch 
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Design  II  Schematic 
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ratio   1%)  will  be  only  slightly  different  for  any  of  the 
four  mechanisms  for  small  values  of    i/D.   The  greatest 
difference  will  occur  for  large    vD.   Table  A-l  shows 
the  value  of  Xc/R  for  a   VD  =  1.5,  which  is  as  large  a 
pitch  ratio  usually  encountered. 

Design  I  Design  II      Pig.  9-A(a)   Pig  9-A(b) 

.478  .430  .466        .431 

Table  A-l 
Values  of  the  ratio  ~  for  a  pitch  ratio  </$=   1.5 

Since  the  mechanism  of  Figure  9-A(a)  yields  the  simplest 
relationship  between  pitch  and  control  rod  displacement,  it 
is  used  throughout  the  study.   The  overall  results,  however, 
will  apply  with  reasonable  accuracy  to  any  of  the  hub 
mechanisms  shown. 

;„  _£_   =    J&  (39) 

TTD  R 


Por  Design  I, 

■r 
c  T~  c 


Kpx  =  ±_AR  +  Pc  bi     +  tc  ll  (40) 


R  AR 


where 

x  =  pump  control  displacement,  in. 

Kp  =  pump  constant,  in^/in 

Kpx  =  flow  from  the  pump,  in^ 

xc  =  control  rod  displacement,  in 

xc  =  velocity  of  control  rod,  in/sec 

Aft  =  control  rod  piston  area,  in2 
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FIGURE  A-9 
Two  Additional  Hub  Mechanisms 
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FIGURE  A-10 
Design  I  Block  Diagram 
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±CAR  =  piston  flow,  ln'/sec 

Fc   =  control  rod  force,  lb. 

KL   =  leakage  constant  of  pump,  in5/!6 

Pc  ^c=  leakage  flow 

Fc   =  time  derivative  of  control  rod  force,  lb/sec 
Kc   =  —  =  compliance  constant  of  the  system  under 


B   pressure,  in5/ib 


V    =  volume  of  oil  under  pressure,  in^ 

B    =  ^        =   bulk  modulus  of  hydralic  fluid,  in 
P  lb 

P    =  hydralic  pressure,  lb/in^ 

J*  Kc  =  compliance  flow,  in^/sec         • 


Equation  (40)  neglects  leakage  around  the  control  rod 
piston  and  expansion  of  the  hydralic  piping.   These  factors 
can  easily  be  accounted  for  increasing  K-r  to  account  for 
piston  leakage3  and  increasing  K_  to  account  for  piping 
expansion.   The  compression  of  entrained  air  bubbles  in  the 
hydralic  fluid  can  also  be  accounted  for  by  increasing  Kc 
appropriately. 

For  Design  I, 

where  Fc   is  that  component  of  control  rod  force  transmitted 
from  spindle  torque,  lb. 

Using  equations  (39),  (40),  and  (41)  the  expression 
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for  the  dynamic  motion  of  the  pump -propeller  system  is 

Input     Ram  Flow  Leakage  Plow 

Kp#  =  ^J^|  A*   4-  %=■  r*-,?^    <eV^_  K^ 

c^putce  a.*»       ,  ^ , . *       A*  cat  bnr- 

^S  J 

(42) 

where  R  =  crank  radius,   ft. 

f^  =  control  rod  ram  sliding  friction,  lb 

/v^  =  control  rod  and  piston,  viscous  friction 
coefficient,  lb-sec 
ft 

Mc  =  control  rod  effective  mass,  lb-sec^ 

ft 


For  Design  II,  as  shown  in  Figure  A-8  the  equation  for 
flow  to  the  power  piston  is: 

:x  "  *cAR  +  pc- 


(y-xc)  L  Kx  -  ±PAP  +  F05L  +  Fn  *S.  (43) 


AR 


Ar 


Where 

y  =  pilot  piston  displacement,  in 
L  =  inlet  port  width 
K_  =  valve  constant,  ln^/in 

The  control  rod  force  in  this  case  is  the  same  as  with 
Design  I  (the  values  of  parameters  will  vary)  with  the  exception 

-so- 


of  an  additional  term  contributed  by  the  safety  spring. 
The  resulting  equation  for  Design  II  is: 

Mechanical 
Input  Ram  Flow  Feedback  Leakage  Flow 

R  — R         w     "  ^P-        £(?.        J 

SPr.^  +cwi         CoVApW^te    -flaw 

^     /  «r~  kl  ^     Spy  iV»4  +«iriv 


Where 
zK  =  spring  constant,  lb/in 


^R  =  spring  force  when  a  =  0,  lb 


C. 


(44) 


The  pilot  valve-stroke  piston  system  shown  in  Figure 
A-7  for  Design  I  is  standard  and  is  fully  analyzed  in 

Reference  33.   It  can  easily  be  shown  that  the  time  constant 

p 
of  such  a  system,  loaded  with  0.4  lb-sec   mass  and  supplied 

with  1000  lb/in  pressure  oil  will  have  time  constants  less 

than  0.1  sec  even  with  excessive  volumes  of  oil  under 
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pressure.   Since  the  transients  of  this  portion  of  the 
system  decay  in  a  small  fraction  of  the  remainder  of  the 
system  transients,  the  pilpt-valve-stroke  piston  system 
can  be  viewed  as  a  simple  gain. 

The  feedback  provided  by  the  mechanical  linkages  of 
Design  I  are  shown  in  Figure  A-10. 

Prom  Figure  A-10  it  is  evident  that 

x  =  &-  (  y  SL    -  xce      \ 


(45) 


Substituting  equation  (45)  into  (42) \   using  relation 


(59), 

Input 


Ram  Flow 


Mechanical  Feedback 


4-    Ku 

7^7 


t  iz.fr  Afr.       >2  Re  kp  ^-l    /p 


(fr-*^m*,    frM"3    _    0<»*O      L 


-/pu 


R 


-  00 


+ 


M,S) 


_  ^W  _  (1^+^^ 


"+  /f  ■  ^-ii^ 


4-    &.  /«M+/J   £) 


/fr-R. 


(46) 
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With  the  exception  of  the  extra  terms  envolving 
the  spring  force  in  equation  (44)  equations  (46)  and 
(44)  are  of  exactly  the  same  form.   They  are  arranged 
so  that  the  analogue  can  be  easily  drawn.   The  analogue 
diagram  for  Design  I  is  obtained  from  the  analogue 
diagram  of  Design  II  by  simply  removing  the  two  spring 
force  feedback  loops  and  adjusting  parameters  to  match 
the  actual  system. 
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APPENDIX  B 
Determination  of  System  Parameters 

1.   Ship  Characteristics 

The  Victory  Ship  (AP2)  was  chosen  as  a  prototype  to 
obtain  representative  parameters  for  the  system.   The  Troost 
B4-55  propeller  curves  were  used  to  determine  the  propeller 
operating  point  for  the  "Victory  Ship  specifications.   Values 
for  wake  fraction  were  taken  from  references  16  and  17. 
These  values  are  summarized  as  follows: 
Ship,  victory,  AP2 

SHP  6000  hp . 

D,  propeller  diameter  18.25  ft. 

Speed,  V 15*5  knots 

RPM,  n  x  60  100  rpm 

Number  of  propeller  blades,  N 4 

Effective  horsepower,  EHP 4650  hp. 

Thrust,  To  9b, 000  lb. 

Propeller  Torque,   Qp , 3.1xl05  lb-ft 

Thrust  wake  fraction,  ujtt     0.348 

Torque  wake  fraction, uugj 0.368 

Thrust  coefficient,  KT 0.310 

Torque  coefficient ,  KQ  0.054 

Pitch  ratio,  f/D  1.15 


Mean  width  ratio,  MWR 0.259 

Speed  of  advance,  Vo  . ... 16.7  ft/sec 

2.   Propeller  Characteristics 

The  Troost  curves  for  Kq  were  then  linearized  about  the 
operating  point.   The  resulting  linear  approximation  to  these 
curves  are  shown  in  Figure  B-l  and  B-2. 

Values  for  K^,  K2>  K,,  L,  K,-,  and  Kg  were  then  obtained 
by  substitution  of  the  appropriate  values  into  equations  (7)> 
(8),  (9),  (14),  (15),  and  (16)  respectively. 

&L  (1.02A)(l8t24)*  lb.seo2  =  288  lb.se02  (7, 

1   4"  bQ         4o7T^ 


K, 


=  ^Qpp4  =  (.476) (1.024) (18. 24)4  lb_sec2 


2TTbQ  201 


.2 


=  863  lb-sec*  (8) 


K3  =  hlUL-     =   .3^1.024)(18.25)5   lb.sec2.ft 
4TT2bQ  4°TT2 

=  1840  lb-sec2=ft  (9) 


K     =     PS?  _    (1.024)  (18.25)-*        lb-sec2 

4TT2b  Tl2   <    -2>  « 

=  69.5  lb-sec2-ft  (14) 
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K5  =     "^pP3         =    (.975)(1.024)(l8.25)3       lb-sec2 
2TTbT  2  7T   (2.5)  ft 

=  420  lb-se°2  (15) 

ft 


k6  =   iix^I   =  (o>iMi.Q24Hie,2gr     lb.sec2 

6    4TT2bT        9.2  TT* 

=  125.2  lb-sec2  (16) 

Curves  of  Kg  were  not  available  for  any  propeller. 
Rupp's  data  for  control  rod  forces  on  a  controllable  pitch 

propeller  while  operating  on  a  tug  (YTB-502)  In  free  route 

4 
were  analyzed  as  follows: 

An  approximate  wake  fraction  was  computed  from  a  formula 

relating  the  block  coefficient  of  the  tug  to  the  wake 

5 
fraction.    The  propeller  speed  of  advance,  Vo,  was  then 

computed  for  the  range  of  speeds  tested.   The  speed  of 


4.  Lewis  A.  Rupp,  "Controllable  Pitch  Propellers," 
Transations  of  the  Society  of  Naval  Architects  and 
Marine  Engineers,  Vol.  56,  1948.  pp.  272. 

5.  S.  Curtis  Powell  and  W.V.  Bassett,  "Practical  Aspects 
of  Torsional  Vibration  in  Marine  Geared-Turbine 
Propulsion  Units,"  Society  of  Naval  Architects  and 
Marine  Engineers,  March  1944. 
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advance,  J,  was  then  computed  from  its  definition  and  the 
available  data  for  each  pitch  ratio  and  speed.   The  corres- 
ponding control  rod  force  was  taken  from  curves  given  in 
Rupp's  paper.   Since  these  forces  represent  both  the 
hydrodynamlc  and  oentrifugal  components  of  torque  converted 
to  force,  the  centrifugal  force  component  of  torque  was 
subtracted  from  the  total  torque  obtained  from  this  force. 
The  blade  torque  coefficient  for  each  pitch  ratio  and 
advance  coefficient  was  then  computed  from  its  definition, 
equation  (19).   The  results  were  then  plotted  and  found  to 
be  very  nearly  linear  as  shown  in  Figure  B-3.   The  tabula- 
tion of  the  calculations  leading  to  the  curves  of  Kg  are 
shown  in  Table  B-l. 

As  stated  in  Appendix  A,  the  blade  hydrodynamic  spindle 
torque  is  a  strong  function  of  spindle  position.   The  actual 
values  for  Kg  could  be  shifted  up  and  down  the  family  of  Kg 
curves  by  shifting  this  spindle  axis.   For  the  purpose  of 
this  study  the  Kg  family  of  curves  were  assumed  to  be 
oriented  such  that  the  curve  for  zero  torque  intersected 
the  origin  of  the  pitch  ratio  -  advance  coefficient  plot. 
This  can  be  Justified  only  in  that  a  propeller  could  be 
made  to  give  this  result.   This  design  would  give  a  much 
more  reasonable  torque  In  the  region  of  normal  operation 
than  merely  using  the  tug  propeller  curves.   The  tug 
propeller  was  designed  for  a  pitch  ratio  of  about  0.7  in 
lieu  of  the  Victory  Ship  with  a  f/D  ±     1.15.   The  slope 
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FIGURE  B-3 
Derived  Kfi  Curves  for  YTB-502 
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TABLE  B-l 

Tabulation  of  YTB§02  Free  Route  Trial 
Calculations  Leading  to  Kg  Curves 

CB  =  JL.  =  333  x  25 =  0.527  =  Block  coefficient 

LBD    (90.64) (25) (9.72) 


ur=  °B  -  o  10  =  t-*2^  -  0.10  =  0.164  =  wake  fraction 
2     '       2 


RPM 
Vs»knots   V9ft/sec    p=5.8  p=6.7   p=7.6    p=8.5 


5 

7.05 

82 

- 

- 

- 

6 

8.46 

100 

89 

80 

- 

7 

9.89 

118 

103 

92 

84 

8 

11.30 

136 

120 

108 

99 

9 

12.70 

155 

137 

123 

111 

10 

14.10 

176 

155 

140 

126 

11 

15.51 

200 

175 

160 

143 

12 

16.92 

- 

200 

182 

166 

13 

18.35 

J 

-Ift 

200 

5 

.544 

,. 

_ 

_ 

6 

.535 

.600 

.668 

- 

7 

.529 

.606 

.679 

.743 

8 

.524 

.595 

.661 

.721 

9 

.516 

.585 

.652 

.723 

10 

.505 

.574 

.635 

.706 

11 

.491 

.560 

.614 

.687 

12 

- 

.536 

.588 

.644 

13 

A 

iL 

^Fc 

.581 

////#/ 

FIGURE  B-5 

YTB-502  Hub  Mechanl 
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P/D 


5 
6 
7 
8 

9 
10 
11 
12 
13 


.610 


Pc  z  10-3  lbs  (includes 
centrifugal  force  component) 

.705    .800    .695 


8.8 

. 

. 

9.5 

9.0 

9.0 

- 

11.2 

9.9 

9.4 

9.1 

13.6 

11.1 

9.9 

9.5 

16.5 

12.7 

10.6 

9.6 

20.1 

14.4 

11.6 

10.1 

22.5 

16.8 

13.1 

11.3 

- 

20.8 

15.5 

13.2 

- 

- 

- 

17.5 

QBH  ft-lb-xlO'^,  centrifugal 
force  component  subtracted 


5 

6 

7 

8 

9 

10 

11 

12 

13 


p/l 


5 

6 

7 

8 

9 
10 
11 
12 
13 


4.15 

- 

- 

- 

4.38 

4.20 

4.24 

- 

5.06 

4.50 

4.33 

4.22 

6.05 

4.91 

4.46 

4.28 

7.28 

5.51 

4.54 

4.16 

8.75 

6.37 

4.82 

4.20 

9.20 

7.31 

5.25 

4.47 

- 

9.00 

6.00 

5.14 

- 

- 

- 

6.51 

K-n 

_  Qb 

Note: p =1.024 
to  match 

AB 

Pn2D5 

Troost's 

curves 

.610 

.705 

.800 

t8?5 

.0278 

_ 

.01915 

.0238 

.0302 

- 

.01604 

.01922 

.0233 

.0272 

.01481 

.01550 

.01740 

.0198 

.01379 

.01340 

.01354 

.01538 

.01291 

.01207 

.01112 

.01200 

.00938 

.01083 

.00930 

.00994 

— 

.00910 

.00820 

.00845 

- 

- 

- 

.00665 
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(m    )  and  spacing  (br>)  were  assumed  to  remain  the  same. 
B 

The  resulting  linearized  approximation  is  shown  in  Figure 
B-4. 

Values  for  K7,  K3,  and  KQ  were  then  obtained  by 
substitution  of  the  appropriate  constants  into  equations 
(22),  (23),  and  (24)  as  follows: 

K    44"  =  tlf&lti^*.  M  Usrf  (22) 

7   4]T2bB     4TT<6)  ft 


K8  =  "»f  **   =  (.650)d.024)(;8.25)4=  ^  ^^    ^ 
2TTbB  2TT(6) 


Kq  =  hBPP5   =  (Q)(1.024)(l8.2g)5  =  Q  (g4J 

47T*bB         4JT2  (6) 


Sliding  friction,  as  represented  by  equation  (25), 
contains  the  thrust  perpendicular  to  the  spindle  axis.  This 
is  the  resultant  of  force  vectors  in  the  plane  of  the 
propeller  disc  as  well  as  perpendicular  to  it.   These 
vectors  change  with  pitch,  speed  of  advance,  and  RPM.  These 
dependencies  are  difficult  to  account  for,  and  Tl  was 
assumed  equal  to  "To  and  constant.   The  dimensions  assumed 
for  bearing  length  and  diameter  are  based  on  dimensions 
scaled  from  existing  CCP  hub  mechanisms  as  shown  in  Ref .  11. 
In  each  case  the  dimension  was  assumed  to  be  a  value  giving 
the  least  final  value  for  sliding  friction  force.   A  friction 
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coefficient,  f ,  of  0.10,  was  assumed. 

QBS  =  fr  TlJ^H    =  (0.10)(l.l4)(9.8xl04)(4.56)4       (25) 
b  (1.4) 


QBS  *  L455X106  lb-ft. 

Since  Q^s  is  so  large,  even  with  low  estimates  for  f , 
r,  l,and  \±    and  a  high  estimate  b,  the  value  for  sliding 
friction  contributed  by  the  control  rod  and  piston  was 
assumed  zero.  i.e.  f R  =  0 

A  sketch  of  the  hub  mechanism  with  dimensions  is  shown 
in  Figure  B-6. 

Reference  11  contains  a  plot  of  control  rod  force  for 
the  transient  conditions  of  going  from  full  ahead  to  full 
astern  while  the  propeller  was  spinning  in  air.  This  curve 
is  shown  in  Figure  B-7.  For  a  spin  test  in  air  the  thrust 
forces  causing  high  values  of  sliding  friction  are  small. 
The  long  region  of  both  transient  curves  that  are  parallel 
with  the  steady  state  (centrifugal  component  only)  corres- 
pond to  the  condition  of  maximum  pump  output  and  constant 
speed  of  the  piston.   The  deviation  from  the  steady  state 
curve  in  both  cases  is  approximately  3  x  104  lbs. 
Reference  U  states  that  the  total  control  rod  travel  is 
0.318  ft.  and  that  complete  reversal  of  pitch  takes  five 
seconds  when  the  propeller  is  spinning  in  air. 

Assuming  that  this  entire  force  is  caused  by  viscous 
friction,  and  the  maximum  ram  speed  is  equal  to  the 
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FIGURE  B-6 
Propeller  Mechanisms  Showing  Only  One  Blade 
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average  speed, 


Xc  =  *2ig,  =  .0636  ft/sec 

^4lb 


the  viscous  force  a  3.0  x  lCT-10  =  ^t/  £( 

)-S( 

ft 


^U  =  4.71  x  105  : 


where  jj  =  total  viscous  damping  coefficient, 


ft 

Assuming  the  viscous  forces  scale  as  the  area  and 
using  the  viscous  friction  term  of  equation  (42), 

9L  (^B*   +>4R,    =    /DvictorjrN   2  R        (4?) 
TTD    -f"       *  '     X^(D  YTB-502/ 

^T   f££  +  y^H)  =  4.71  x  105   (I2*ffi  >2  1#2 
1111     R  9.5 


i-  (<B^  +  /^RR)  =  3.65  x  105  lb-sec 
TTD    r  ft 


In  order  to  obtain  a  value  for  the  polar  mass 

moment  of  inertia  of  the  blade  about  its  spindle  axis,  the 

blade  is  assumed  to  be  an  elliptical  plate  with  an  average 

thickness  ss  ft. 
6 

The  formula  for  the  moment  of  inertia  of  an  ellipse 
about  its  major  axis: 

JiBb  "  %     £^J      m  (48) 

8T21 
where  t  =  thickness  of  plate,  ft. 
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(%         =  density  of  metal,  lb-sec2 
r  ft^-ft 

a     =  length  of  minor  axis,  ft, 

b     =  length  of  major  axis,  ft. 

**iBb  =  mass  moment  ot   inertia  of  blade  about  major 
axis,  lb-ft-sec2 

The  axis  of  the  spindle  is  offset  so  that  approxi- 
mately 40$  of  the  area  is  forward  of  the  spindle  axis. 

/.  JiBH  =  tfe  a3/2blT  N  +02tpTTabN 


where 

C  =  distance  of  spindle  axis  from  the  major  axis,  ft. 
N  =  number  of  blades 
Substituting  values : 

jiBu  =  J  (522  ,  (g.y»(9-i2S)  *  (4Q) 

*         6  32. 2)  8T21  (  9' 

+  (.67)2  I  &H   )  TT  (5.37)  (9.125)4 

6   32.2 

JiBN  =  808.5  lb-ft-sec2 

Estimating  the  mass   of  the   control  rod   a  length  of  20 
feet  was  assumed. 

Mc  =    PcTT(Rc   )2  L  =  it22      (li6S   ^2    (20)   =  740  lb-sec2 
*  12  32.2        12      ;  ft 

where 

o 

Pc  =  mass  density  of  steel,  lb-sece 

R„  =  radius  of  control  rod,  7.67  inches 


L  =  rod  length,  ft. 

This  is  a  large  estimate  for  control  rod  mass,  but  it 
is  shown  in  Appendix  C  that  inertia  loading  can  be  neglected 
even  for  this  large  estimate. 

Allowing  a  30%   increase  for  entrained  water  an  estimate 

'  2 

of  JiB  =  1050.0  lb-ftsec  is  obtained. 

The  assumption  of  an  elliptical  plate  was  also  made  in 
evaluating  the  constants  for  centrifugal  blade  torque,  K10 
and  K11# 

The  formula  for  an  ellipse: 

7?       y2  _  1 

a2  +  b2  "  ?  '  (50) 

Translating  the  axis  to  the  blade  root  and  solving  for  X: 


X  =  a  (1  -  if  )*  +  k 


%    "fe  >  +f      •  <51) 

The  blade  section  area,  A,  is  a  function  of  distance  from 
the  root, 


A  =  t2X  =  2t 


U   b2  }   +  2 


(52) 


The  distance  from  the  center  of  the  area  to  the  spindle 
axis  in  the  x  direction  is  assumed  constant  =  c.   The 
distance  of  the  center  of  the  blade  section  to  the  spindle 
axis  in  the  y  direction  is  a  function  of  rake. 

y  =  mr^  (53) 

where  nip  =  slope   of  blade  from  vertical. 
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Substituting  equations  (52)  and  (53)  into  equation  (36): 
tip 


K10  " 


»f_     \      AxydJ? 


(36) 


K10  .  I  P 


hub 

p  ,  *i?  r  a(i  -£)*♦* 


2tL 

"J  hub 


b2 


5mr 


M 


If  nip  is  small  the  vertical  distance  is  approximately 

equal  to  the  blade  length, 

Atip 
K10  =  ■  PP  \  2t  fa  (J  -  \    )*  +  |  J  cn^IdY  . 


(54) 


hub 


Integrating: 


K10  =  N  PP  ^tci^a 


+  N  p  4tefflr  b 


I  (b2 


v 

2 


I2)*  ♦  g-  Sin  "til  )  |  2 


b 
2 


Assuming  nip  =  i  and  substituting  values  for  the 
remaining  variables: 


K10  =  *  (55*  )  (4  )  (.67)(i)(2.685)   2^12^" 


32.2    6 


(2) (B)(4) 


+  4  (552.  )  (i)(.67)(i) 


32.2 


(9.125)3 
(64) 


Kin  =  105.7     lb-ft-sec2 


L10 
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The  formula  for  the  moment  of  inertia  of  a  rectangle 
about  it  centerline: 


UMIF 


12 


where  the  length 
also 


=  X 


LMAX 


(X)3t 
12 


(55) 


(56) 


Using  formula  (51)  for  X,  the  expression  for  K-^  can 
be  approximated  using  equation  (37). 

ptip 

K11  =  NnDfp    K  (i^  -  y  di. 

^root 


(37) 


K 


11 


nttd  P  . 


ife 


b2    )      + 


H 


12 


12 


i5  [.(*-S>**i] 

12     L  J 


dY 


(57) 


After   expanding  the   cubed   term,    (57)    integrates   to 


K 


11 


=  NTTDfLtab 


12 


5?     +  8 


3b£~| 
32  J 


NlTD^    t5a 


£D"7     -  NTTDPt^b2 
32    I  IE. 


24 


(58) 


Expression    (58)   reduces   to 
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Kll  =  Hof>f  bt    (   a2  r2+T(6b2   n      ^1         (59) 

768a    i  L  J       J 

Upon  substitution  of  the  limits. 

Substituting  values  for  the  various  parameters: 

K„  =  4TT  (18.25)  CJ£g)(9.125)<fc)  (&.*)*[&   37)2 
768  (5.37)         C      L 

+  TT6(9.125)2  -  3lf]-32   (9.125)(^)2> 
Kll  =  945  lD"sec2 


3.  Pitch  Control  Mechanism  Parameters. 

The  value  for  R  is  obtained  from  the  geometry  of  the 
hub  mechanism  and  is  shown  in  Figure  B-6  to  be  about  1.2  ft., 
which  is  the  maximum  feasible  for  this  size  propeller. 

A  maximum  operating  pressure  of  1000  (b/in^  was  assumed. 

The  piston  area  required  is  obtained  from  summing  the 
forces  on  the  control  rod  and  dividing  by  the  pressure. 

Combining  equations  (21)  and  (35)  and  substituting 
operating  point  valves  for  the  variables  yields  the  spindle 
torque  cause  by  hydrodynamic  and  centrifugal  forces. 

Qbh  +  Qbc  =  <K7  -  Kn)  p^2  -  ¥0^  "  <K9  +  Kio)aJ  2 

QBH  +  QBC  =    (482~  945) (21.0)    (ipj  )2 

-  1955   (16.7)  (±22)   -    (0  +  105.7)  (i§I) 

3  3 

QBH  +  QBC  =  1*063  x  107  lb-ft 
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when  the  maximum  torque  condition  prevails. 

QBS  =  +  frTx^  g   -  (.10)(l.l4)(9.8xl04)(4.56)4 
b  1.4 


QBS  =  -     1.455  x  106  lb-ft 

If  these  were  the  only  torques  to  overcome,  le 
neglecting  viscous  friction  and  inertia  loading,  the  total 
torque  would  be  the  sum  of  the  torques  above. 

Qbh  +  Qbc  +  Qbs  ■  1'063  x  lo7  +  1-455  x  lo6 

^BH  +  QB0  +  QBS  =  12.09  x  106  lb-ft 
Dividing  by  R  yields  control  rod  force. 

QBH  +  QBQ  ^  QBS  =  12.09  X  106     lb  .  1#005  x  107  lb 
R  1.2 

The  ram  area,  Ajj,  becomes 

1.005  x  107  lb-in2   =  x  Q  x  1Q4  in2 
1000  lb 

A  piston  of  this  size  would  not  fit  within  the  hub  of 
4  ft.  diameter  for  the  Design  II  type  mechanism  and  is 

excessively  large  in  any  case.   Hence  a  piston  area  of  775 

2 
in  is  assumed  which  gives  a  piston  radius  of  15.75  in.  This 

size  is  about  the  maximum  that  can  be  housed  within  the  hub 

and  for  either  design  there  is  a  likelihood  of  stalling 

under  the  heavier  load  conditions.   A  larger  diameter 

causes  extreme  demands  on  pump  output  in  order  to  obtain  a 
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reasonable  response  for  the  pitch  setting  mechanism.   This 
pump  must  be  very  large  even  with  this  size  piston. 

As  is  shown  in  Appendix  E,  the  time  constant  of  the 
pitch  control  mechanism  is 

Tp  =  RAR 

-4r-  (60) 

KpTTDM 

It  is  necessary  to  reduce  this  time  constant  to  as  small 
a  value  as  possible.  The  larger  this  time  constant,  the  less 
effective  the  control  mechanism  becomes  in  RPM  regulation. 
This  is  shown  graphically  in  Figure  E-6.  As  Tp  increases  the 
frequency  at  which  the  ratio  of  RPM  fluctuation  to  wake 
disturbance  begins  to  increase  for  all  systems  with  regu- 
lation decreases. 

Higher  values  of  Kp  result  in  lower  values  of  Tp.   The 
higher  values  of  Kp  correspond  to  large  pump  capacities. 
Assuming  a  300  H.P.  variable  displacement  is  the  maximum 
available  Kp  is  calculated  as  follows: 

A  300  HP  variable  displacement  pump  is  approximately 
capable  of  a  maximum  delivery  of  120,000  in^/min  =  2000 
in-y  sec. 

The  maximum  control  arm  displacement  is  arbitrary  and 

depends  on  the  length  of  the  lever.   A  reasonable  value 

for  this  parameter  is  4.368  inches. 

Hence, 

K  _  2000   in3      _  458  In? 
4.368  in-sec         in-sec 
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The  leakage  constant,  KL,  can  be  estimated  by- 
assuming  a  maximum  %   of  full  pump  displacement  at  full 

pressure. 

in3 

KL  =  2Q0Q  ^    (assumed  %)  (61) 

1000  ^_ 
in2 

K  =  o.2  == —    for  10#  leakage 
L      sec-lb 

KL  =  0.02  —     for  1%   leakage 
sec-lb 

The  mechanical  feedback  shown  in  Figures  A-7  and  A-10 
is  determined  by  the  length  of  linkages  e  and  ^  or  the 
length  and  width  of  the  valve  port  in  the  case  of  Design 
II.   Increasing  this  feedback  decreases TL,  but  also 
decreases  the  overall  pitch  mechanism  feedback  gain  as 
shown  in  Figures  E-4  and  E-5.   This  corresponds  to  requir- 
ing a  progressively  longer  travel  of  mechanical  input  at 
y  in  Figure  A-7  or  on  A-8.  The  sensing  devices  whether 
tachometers  or  torquemeters  will  produce  electrical  or 
weak  mechanical  outputs.  A  long  mechanical  travel  at  y 
in  either  case  introduces  additional  time  lags  that  reduce 

system  performance  because  of  inertia  loading  at  a  point 

o 
in  the  circuit  having^low  power  level.  Reasonable  values 

for  the  linkage  lengths  were  therefore  assumed  as  follows: 

e  as  5  inches 

f  =  15  inches 
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a  =  5  inches 

b  =  5  inches 

The  G2  which  corresponds  to  the  gain  of  the  stroke 
piston  system  is  assumed  to  be  one,  i.e.  it  acts  purely 
as  a  force  amplifier. 

M  g  1262e*       =  (12?(1?(^(1.2)    m   0<063 
(e  +  fy  )1JD  (5+15)  TT  (18.25) 

The  time  constant  for  the  stroke  piston  system  is 

smaller  than  any  other  time  constants  in  the  system  by 

several  orders  of  magnitude  and  is  therefore  neglected. 

The  compliance  constant,  K_,  is  determined  from 

c 

equation  (62). 

Kc  =  J  (62) 

where 

V  =  volume  of  oil  under  pressure,  in-' 

1  2 
B  =  bulk  modulus  of  hydralic  fluid,  i=- 

lb 

B  =  0.27  x  106  i&   for  standard  hydralic  oils. 

in^ 

The  maximum  piston  travel  is, 

xc    =  R£  =  (1.2)(21)    u   Q#44  ft# 
MAX   TO     TT  (18.25) 

The  volume  of  oil  under  compression  in  the  piston  is, 


6.   Truxal,  J.G.,  op.clt.  t  p.  15-44. 
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AR  xc         =   0.44  x  775  x  12  =  8,200  In3 
*  CMAX 

Assuming  a  10  foot  length  of  piping  under  compression 

with  a  diameter  of  3.4  inches,  the  volume  under  compression 

in  piping  is, 

(20)(12)(^A4)2Tr  =  1180  in3 
4 

The  total  volume  under  compression: 

8,200  +  1,180  =  9,380  in5 

,\  kc  =  |  =  228p_    =  .0336  i£ 
C   B   0.27  x  10°         16 


4.   Power  Plant  Characteristics 

Propeller  inertia  (including  entrained  water), 
propeller  damping,  power  plant  inertia  referred  to  the 
propeller  shaft,  and  engine  damping  assuming  the  turbine 
is  a  constant  torque  device,  is  taken  from  Sutton's 
vibration  analysis  of  the  Victory  Ship,  (AP2). 

J±p  =  3.09  x  104  lb-ft-sec2 


JiE  =  3.31  x  105  lb-ft-sec2 


oCE  =  3.08  x  104  ft-lb-sec 


,4 

XT 

rad 


7.   A.D.  Sutton,  "Torsional  Vibration  of  a  Geared  Turbine 
Propulsion  System."  New  England  Section  Paper  of  the 
Society  of  Naval  Architects  and  Marine  Engineers, 
October  1953,  p. 27. 
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°CV  =   7.54  x  10  ft-lb-scc 

rad 

oCp  is  also  found  in  Appendix  E  to  be  equal  to  K-, 
and  can  also  be  determined  from  propeller  curves. 
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APPENDIX  C 
Development  of  the  Non-Linear  Analogue 

The  final  expressions  derived  in  Appendix  A  could 
be  programed  directly  for  analogue  computation.   Simpli- 
fication of  equation  (46)  greatly  reduces  the  complexity 
of  the  program  and  is  therefore  desirable. 

The  compliance  term  in  equation  (46)  contains 

d_  f  (K7-Kn)  v(Jj2     -  KSVpOJ   -  (K£j!*10iQj2 
dt   /    E  ~B~        H 


1   (f„  ^rTJ  N; 


(63) 


bR 
Performing  the  derivative, 

i  f(K7-K11)(paJ2    +  2CJpOJ  )    -  K8V0^>   -  KqU>V( 

L  +     i_     (f„   +  frTx  X   g      ) 

dt  *  bR 

Cc)  is  the  variable  being  regulated,  and  even  without 
regulation^  is  small.   Therefore,  those  terms  containing 
u)  are  assumed  zero. 

The  term  (fR  +  fpT^X  g   )  ls  pictured  as  a  function 

b 

of  p  in  Figure  C-l. 
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FIGURE  0-1 
Sliding  Friction  as  a  Function  of  p 
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This  sliding  friction  term  is  discontinuous  at  $  =  0. 
The  time  derivative  of  the  friction  term  is  infinite  as  the 
direction  of  motion  of  pitch  is  reversed  and  therefore 
causes  an  Impulse  of  flow  or  a  shock.   It  does  not  allow 
a  volume  to  flow  and  therefore  can  be  neglected  in  the 
analogue.   The  presence  of  the  term  indicates  that  a 
hydralic  shock  will  be  generated  for  each  pitch  reversal 
in  the  actual  system. 

Expression  (63)  reduces  to 

5  [  <*7  -  Kll)  $"2  '  ^0^]  <64> 

with  the  above  simplications. 

Expression  (46)  contains  two  terms  in  ff.   They  are 

Ju    (fiSL     +     MCR)   ^L  , 

AR 


and 


L     &L.     +JLR  )  h. 


Substituting  values  into  these  two  terms : 


L     \222*1     +  740(1.2) 


.0336  b' 

-~    =_r—       (52.5) 


TTD         1.2  J       AR  TTlAr 

and 

|_     (2.08  x  10?)   £L     =  I (4.16  x  107) 

TTD  ar       ttdar 
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Hence  the  first  term,  envolving  inertia  and  leakage 
can  be  neglected. 

Substituting  expression  (47)  for  expression  (46)  in 
equation  (37)  results  in  the  following  three  terms 
containing  p. 

J  J       (1,2) (12) (775)  i   =  194.9  p  in3/sec 
TTD   R     TT(18.25)  ' 

p  ^ ( °£^  +yRR  )  =  (.2)(1.742  x  IP6)   . 

Fl)Aa       R  (18.25)  (775)   P 

=  0.785  p 
(For  10#  leakage) 

^  Kc(*r*ll>OJ  2  =  (>Q336)  (480-945)  (^J)2 
AR  775  * 

=  2.20  p   in^/sec 


Hence,  ignoring  these  last  two  terms,  equation  (46) 
reduces  to: 

R 


V  =  *  #U  ar  +  —  I  (K7"Kn)  ml    -  K8r0uJ 

AR   L  R      —r  — 


-(K9  +  K10)  <J_     +  (fR  +  frTjJ  ■  -1  (65) 

R         bR         J 
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Equation  (65)  and  the  remaining  equations  for  Design 
I  derived  in  Appendix  A  are  programed  for  analogue  compu- 
tation as  shown  in  Figure  0-2.   Substitution  of  the 
parameters  determined  in  Appendix  B  is  shown  in  Figure  C-3. 

The  analogue  computer  used  will  not  accept  signals  of 
greater  than  50  volts.   Also,  when  the  output  of  any 
amplifier  exceeds  50  volts  the  amplifier  response  becomes 
non-linear  and  saturates  completely  at  about  57  volts. 
For  this  reason  it  was  necessary  to  use  the  following 
scales  for  variables : 

pitch,  1  f t  =  0.1  volts 

shaft  speed,  1  ESi  =  0.4  volts 

sec 

speed  of  advance,  lO—  =0.1  volts 

sec 

A  scale  larger  than  this  saturates  the  path  containing 

2 
p  co     terms . 

Care  must  be  taken  when  using  different  scales  for 
variables  when  non  linear  elements  such  as  the  multipliers 
are  used. 

The  general  rule  for  multipliers  is  to  divide  the 
amplification  following  the  multiplier  by  the  product  of 
the  two  scale  factors  of  the  input  variables. 

Scaling  through  linear  elements  such  as  amplifiers 
and  integrators  is  accomplished  by  increasing  the  coefficient 
for  each  input  path  containing  pitch  and  speed  of  advance 
as  a  variable  by  a  factor  of  four. 
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Finally  It  Is  necessary  to  move  various  gains  around 
the  system,  following  the  rules  of  block  diagram  trans- 
formation, so  that  individual  amplifiers  will  not 
saturate  and  the  inputs  to  each  amplifier  do  not  exceed 
four,  all  with  coefficients  having  their  decimal  point 
within  three  digits  of  each  other.  An  example  of  this 
last  step  is  shown  in  Figure  C-4-.  The  gains  must  also 
be  arranged  so  that  no  amplifier  is  called  upon  to 
provide  a  gain  greater  than  10  • 

The  portion  of  the  analogue  marked  "Selector- 
Amplifier  Section"  effectively  produces  the  response 
pictured  in  Figure  C-l.  The  bias  and  amplifier  gains 
are  set  for  a  very  small  sinusoidal  input  to  give  the 

correct  voltages  corresponding  to  the  term,/fR  +  frTxx  N  j 

bR   /  ' 

Figure  0-6  is  the  analogue  with  all  scaling 
completed  and  the  program  ready  for  the  computer. 
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FIGURE  0-4 
Example  of  a  Block  Diagram  Transformation 
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APPENDIX  D 
Non-Linear  Analogue  Computation  Results 

1.  Difficulties. 

The  Philbrick  Analogue  Computer  of  the  Mechanical 
Department  Systems  Project  Laboratory  at  M.I.T.  was  used 
to  obtain  data.   This  computer  utilizes  electronic 
multipliers  to  perform  multiplication  and  division.   The 
multipliers  produced  a  hum  that  was  within  the  band  pass 

of  the  amplifiers.  Although  correct  average  values  for 

2        •  2 

the  variables,  60  »  V0o>,  V0<^,  and  p6c>  were  obtained,  the 

hum  superimposed  on  the  signal  was  transmitted  through 

each  amplifier  and  in  some  cases  obscured  the  data.   The 

power  plant-propeller  system  and  the  pitch  control 

mechanism  were  operated  separately  with  satisfactory  results, 

but  the  connecting  link,  marked  G  in  Figure  C-2  introduced 

hum  for  the  variable  p  saturating  the  integrator  that 

followed.   For  this  reason,  no  overall  system  response  was 

obtained. 

2.  Useful  Results. 

The  response  of  the  power  plant-propeller  system  to  a 
sinusoidal  input  of  VQ  was  obtained.   A  two  channel  Sanborn 
recorder  was  used  with  the  input,  V0,  on  one  channel  and  the 
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Ca)     response  on  the  other.   The  sSanborn  recorder  stylus 
deflection  was  calibrated  against  a  computer  amplifier. 
Each  amplifier  has  a  provision  for  adding  up  to  1  50.0 
volts  to  the  output  in  increments  of  0.1  volts ♦  With 
one  volt  per  cm.  deflection  the  recorder  output  was 
linear  over  t   2  cm.   The  maximum  output  recorded  was 
1  2.2  cm.  and  the  accuracy  of  readings  is  approximately 
0.1  cm.  The  results  are  shown  in  Table  D-l  and  are  plotted 
in  Figure  E-6. 

The  time  constant  of  the  power  plant-propeller 
system  was  observed  to  be  approximately  3.28  sec.   This  was 
determined  by  injecting  a  step  of  Q  and  observing  the 
response  in  6J  .   The  decay  to  steady  state  in  to    occured 
in  about  10  sec.  or  three  time  constants. 

The  values  of  QB  and  Qp  generated  had  an  average 
value  within  5%   (3.1  volts  and  7.0  volts  respectively)  of 
the  calculated  values  at  the  operating  point,  but  a  sinu- 
soidal hum  of  volts  was  superimposed  on  both. 

The  node  corresponding  to  the  variable  p*  had  a  sinu- 
soidal amplitude  of  .6  volts  for  a  8  volt  sinusoidal  input 
for  V0  with  no  friction  feedback.   With  friction  added  this 
signal  was  reduced  to  less  0.1  volts  (the  approximate  level 
of  noise).   This  is  an  indication  of  the  importance  of  the 
strong  stabilizing  effect  of  sliding  friction. 

The  feedforward  path  containing  v"0  cO    had  no  effect  on 
p  even  with  an  input  of  30  volts  for  AV0. 
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TABLE  D-l 
Power  Plant-Propeller  System  Frequency  Response 


frequency 
cycles/sec 


0 

.01 
.02 
.03 
.04 
.05 
.06 
.07 
.08 
.09 
.1 
.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 
1.0 
2.0 
3.0 
4.0 
5.0 


AV0 

Loo 

A                         \ 

Input 

Output 
(Volt  Sxi) 

AUJ 

101  o 

g       \bU) 

(Volts) 

AV0    1 

idpr; 

3.0 

0.15 

0.050 

• 

13.0 

3.0 

0.15 

0.050 

- 

13.0 

3.0 

0.15 

0.050 

- 

13.0 

3.0 

0.20 

0.066 

• 

11.8 

3.0 

0.20 

0.066 

- 

11.8 

3.0 

0.20 

0.066 

- 

11.8 

3.0 

0.20 

0.066 

- 

11.8 

3.0 

0.20 

0.066 

- 

11.8 

3.0 

0.20 

0.066 

- 

11.8 

4.0 

0.20 

0.050 

- 

13.0 

4.0 

0.20 

0.050 

- 

13.0 

8.0 

0.38 

0.048 

- 

13.2 

11.0 

0.62 

0.056 

- 

13.4 

15.0 

0.84 

0.056 

- 

13.4 

18.0 

0.79 

0.044 

- 

13.5 

24.0 

1.03 

0.043 

- 

13.6 

26.0 

1.10 

0.042 

- 

13.7 

28.0 

1.16 

0.0415 

- 

13.8 

30.0 

1.25 

0.0415 

- 

13.8 

38.0 

1.52 

0.0400 

- 

14.0 

38.0 

1.45 

0.0380 

- 

14.2 

38.0 

1.0 

0.0264 

- 

15.8 

38.0 

0.50 

0.0130 

- 

18.8 

40.0 

0.00 

0 

- 

00 
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APPENDIX  E 
Linear  Incremental  Analysis 

1.   Simplification  of  the  model  equation  (46)  contains  the 
term, 


12G2  eRK 
(e  +  fj  )7TD 

which  is  mechanical  feedback. 

Define  12GQeR  ,  g-^ 

2        a  M  (66) 

(e  +  fy  )7TD 

As  shown  in  Appendix  C,  a  reasonable  value  for  this 
term  is 

M  =  0.063 

As  shown  in  Appendix  D,  the  terms  P   K"RN    i  /  ^  v 

7TD  (-|-  +  >HaR)E£.  , 
*         AR 

and  _J^     o    9  are  small  in  comparison  to  the  remaining 
AE    R 

terms  of  equation  (45).   Neglecting  these  terms  the 
analogue  for  the  pitch  setting  mechanism,  as  shown  in 
Appendix  B,  can  be  reduced  to  that  shown  in  Figure  E-l. 
This  is  equivalent  to  neglecting  the  corresponding  term  of 
equation  (46),  which  reduces  to 


#dAr  t 


(K7-Kll}  i^l     -  K8V0 
R     "S 


-(K9  +  Kio)  u£       ;   (fR  +  frT^  N  )  (67) 

R  *>R 
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FIGURE  E-l 
Simplified  Pitch  Control  Analogue 
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FIGURE  E-2 
Pitch  Mechanism  Analogue  Neglecting  Leakage 
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If  f£  is  sufficiently  small  the  leakage  feedback 
AR 

shown  in  Figure  E-l  can  be  neglected  and  the  system  reduces 
to  that  shown  in  Figure  E-2. 

The  transfer  function  for  the  system  shown  in  Figure 
E-2  is 

iisi-   Kp7TD 
yTsT  ~     ^r 


S  +  -* (68) 

RAr 

This  constitutes  a  type  zero  system  having  a  simple 

Q 

lag  with  time  constant. 
T  =  % 

~~^^  (69) 

Substituting  values  determined  in  Appendix  C, 

rp  =   (1.2)  (775) =  .563  sec. 

(45y)TT(l8.25)(.063) 

Leakage  has  the  effect  of  reducing  this  time  constant. 

The  two  largest  contributors  are  friction  and  the  pitch 

component  of  blade  torque.  Figure  E-3  shows  these 

additional  feedback  paths  added  to  the  system  shown  in 

Figure  E-2. 

p 
The  system  transfer  function  that  includes  the  poo 

feedback  term  is 


8.   J.G.  Truxal,  Control  Engineers  Handbook.   The  Maple 
Press  Co.,  York,  Pa.  1958,  p. 2-35. 
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PIGURE  E-3 

Pitch  Mechanism  Analogue  With  Major  Leakage 
Factors  Included 
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KpTTD 

Pis)    =     AfiS P_  (70) 

yTs)  S   +  KpMTTDRAR  +    (K7=Kn  -,  )  TTD^KL  ■  v 

and  now,  * 

A  2R2 

rp=  fsJ s . (7D 


KpMTTDRAR  +    (K?  -  K^TTcO2  K  D 


TP  = 


L4 

(775)2  (1.2)2 

(458)  (,063)  TT  (18.25)  (1.2)  (775 )  + (276-482)^)* (L8.25X2) 

7^  =  0.541  sec. 

The  computation  above  is  based  on  a  leakage  of  10$  of 
pump  output  at  maximum  flow.  The  leakage  caused  by  friction 
can  be  approximated  as  reducing  the  pump  output  or  effectively 
reducing  Kp  as  shown  diagramatically  in  Figure  E-4- 

The  effective  Kp  becomes: 

K.   =  KP*max  -    <VS£=LJ     £L 

KP   br      AR  (?2) 

XMAX 

Ki  _  (458)  (4.368)  -  (1.210  x  106)  &* 
p  (4.368) 


k'  =  384  ln3 

in-sec 

The  new  Tp  becomes 

r  =  5£&  =  (1*2)  (775) (73) 

K^TTDM   (384)TT  (18.25)  (0.63) 
Tp  =  0.654 
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FIGURE  E-4 
Effect  of  Leakage   on  K. 
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FIGURE  E-5 
Regulation  by      '    Sensing 
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The  bad  effect  of  increasing  7p  is  shown  in  the  gain 
plot  of  Figure  E-6.  As  Tl  increases,  the  gain  defined  as 
the  ratio  of  shaft  speed  variation  to  make  fluctuation 
increases  and  at  progressively  lower  frequencies.   It  is 
expected  that  10$  leakage  is  excessive  and  the  "C  computed 
for  1%   will  be  used  in  all  future  calculations.   That  is 

7p  =  0.560 

Viewing  propeller  torque  as  a  linear  function  of  p, 
u>,  and  V0, 

Qp  =  Kjv0  +  K|p  +  K|60  +  Qp0  (74) 

where 


ipO 


=  operating  point  torque,  ft-lb 


KJ 

-  ^Qp 
^  v0 

p,co 

4 

_  >>QP 
^p 

Vq,^ 

K* 

_  ^>QP 

V0»P 

=  °c 


(75) 

(76) 
(7b) 


Values  for  the  above  K*'s  are  obtained  from  propeller 
curves  by  allowing  the  variable  in  the  denominator  of  the 
partial  derivatives  vary  over  t  .1  of  their  operating 
value,  obtaining  the  Kq  at  each  of  the  two  points  and 
solving  for  the  Qp  and  dividing  the  result  by  0.2. 

Calculations  for  Kj,  Kjp,  and  K-z,  using  Troost  s 
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B4-55  propeller  follow: 

1  "  Tv       — r^ —  (79) 

AVo    (.2)nD 
AKn  =  KQ  (for  J  =  1.1J0)  "  Kq  (for  J  =  0.9J0) 

where  JQ  is  the  operating  point  advance  coefficient. 

.'.  K*  =  (-007)  (1.024)  (^)2(18.25)5  =  _6610  lb 

1      (12) (i§) (18.25)  J 

4  =  77  -  *f¥*^  (so) 

*     (.2)  £  x  D 


where  A  KQ  =  KQ  (for  J  =  1.1  (£)q) 


-  Kn(for  E  =  0.9  (P  )  ) 
^     D        D   o 


*  =  (♦02)(1.02j0(;rf(l8.25)4  =  2.96  x  104  lb. 
2       (.2)(1.15) 


*  ..  A3.  =  (  aKq)  p  n2D5 
*3  "  A  co  An  2TT 


K 


In  this  case  p/D  and  VQ  remain  constant 

but  J  =  _o   changes  as  n  =  —  changes  by  -  0.1. 
nD                2TT 

*  _  (.006) (1.024) (i^)2 (18. 25)5  =  7.54  x  104  ft-lb-sec 


3 


(.2)(|2)  2TT  rad. 
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In  summary: 

K,  =  -6610  lb-sec 

K*  =  2.96  x  104  lb 

K*  =  7.54  x  104  ft-lb-sec 
3  rad. 


Using  primes  (  )  to  Indicate  incremental  quantities, 
equation  (74)  is  incrementally, 

Prom  equation  (10), 

-oCE^-    (K*  V0  ♦  YL\  ,'    +  K%(4  )   =  J^^'  (82) 


Qp  =  K?Vl   +  K*P*    +  kIo)  (81) 


(83) 
(84) 


or  using  the  Laplace  transform, 
60  (s)       (oCE  +   oCp)    +  j^   sj=  K*  p'(s)    +  kJVq(s) 

CO  (a)   =       EgPU)    *  Kl   Vo    <s> 
(cCB+flCP^+  4TS 

The  block  diagram  for  which  is  shown  in  Figure  E-3 


2.   Shaft  Speed  Sensing 

Sensing  angular  velocity  and  adding  the  pitch  setting 
mechanism  to  the  diagram  of  Figure  E-3  is  shown  in  Figure 

E-5. 

Define:  oCn  =  °C  +  *S»  =  3.08xl04+7.54xl04=10.62xl04  ft-Jb-sec 

T   *    p  rad. 
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T=  J^     =  ^Q  l   1q5,    =  3.28  sec. 


^E+°S    10'62  X  ^ 

Q   =  gain  for  angular  speed  sensing 
The  transfer  function  between  V0'(s)  and  60  (s)  is: 

°>'t  tf)  s  -5l  (i-Tps) 

V0  (s)     ■■  I ^ (yc>) 

(f  S+l)(fp+l)  A  ^a> 

oC  TM 

The  selection  of  G^  is  usually  based  on  criteria  for 
transient  response  or  stability.  As  G^,  increases  the 
effectiveness  of  regulation  improves  for  lower  frequencies 
but  more  severe  peaking  occurs  near  resonance.  A  compromise 
value  of  G^  giving  a  damping  factor,  J«  0.4  was  chosen  and 
used  as  the  criteria  for  o>  and  Q„  sensing  designs.  A  G^ 
giving  a  damping  factor  of  one  provides  essentially  no 
improvement  of  regulation  overthe  fixed  pitch  system. 

The  damping  factor  is  defined  for  the  transfer  function 
when  the  denominator  is  placed  in  the  form: 

S2  +  2j?  C0ns   +  COn2  (b6) 

where 

T  =  damping  factor 
LA     =  natural  undamped  frequency,  sec"-1- 

Placing  the  denominator  of  equation  (85)  in  the  form 
of  expression  (86),  and  solving  for  Geo  and  ^  n  with  F=0.M- 


yields: 


G 


CD 


(%+  Tp)2^^  =  (3.28+.560)2(3.08xl04)(.063) 
TETp(.64)K|  '  (3.28)(.560)(.64)(2.96zl04) 

=  0.82  ft-aec  (87) 

rad 

O)     =  K£  G    *  =  2.96  x  104  (.82) 

11  °C/LT*  M    (3.08  z  104)  (.560)  (3.28)  (.063) 

=  2.61  seo'1  (88) 


System  response  for  angular  speed  sensing  is  plotted 
in  Figure  E-6. 

3.   Torque  Sensing 

If  a  torquemeter  were  used  in  lieu  of  a  tachometer  for 
sensing,  regulation  of  shaft  torque  might  be  accomplished. 
This  in  turn  would  regulate  shaft  angular  speed. 

Figure  E-7  shows  the  block  diagram  for  the  system  using 
torque  sensing. 

The  V0  to  t*J  transfer  function  is  now 


2*9  -•    l_r       (rpS+D 

vi<a"  — £j — : — 


ct         My 


(89) 


Using  the  same  criteria  of  f=   0.4  in  the  setting  of 
Gq  leads  to  a  quadratic  equation  in  Gq. 
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FIGURE  B-7 
Block  Diagram  for  Regulation  by  Torque  Sensing 
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Solution  for  GQ  follows: 
Prom  expression  (89), 


^£l 5 =2^nf     .  (90) 


■p 


.560  +  3.28  +  7.54xl04(.565)  +  2.96xl04GQ(3.28) 

10.62x10*    >m063 

(.563)(3.28) 

=  2(.4)con 

2.301  +  8.35  x  105GQ=   .80>n 

2.88  +  1.042  x  106Gq  ■  6U  n  (91) 

also, 

L°Ii  =A  +  7.54  x  10*  +  2.96  x  104GqW  1 \   *      (93) 

10.62  x  104  .063  /    (3.26x.560   / 

°n  =   Cl.708  +  4.70  x  105gq)*  (94) 


Combining  expressions    (91)   and    (94), 

1.708   +  4.70  x  105gq=    (2.88  +  1.042  x  106GQ)2 

Rearranging  terms 

«►  1.09  x  1012GQ2+  5.53  x  106GQ+  6.09  =  0  (95) 
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The  positive  root  of  this  equation  Is  approximately 

+  6.78X|0^  ft    =  Go    This  value  of  G0  Is  not  un- 
lb-ft  H 

reasonably  small  when  It  Is  noted  that  the  steady  state 

5 
shaft  torque  =  3.1  x  10  lb-ft. 

Substituting  the  relevant  positive  value  into  equation 

(91)  yields, 

60n  =   (2.88  +  1.042  x  106  x  6.78  x  10"7)   sec"1 
^>n  =  2.17  sec"1 

The  *2  (g)        is  plotted   in  Figure  E-6  for  these  values 
V0(s) 

of  Gq  and  to    . 


4.   Thrust  Sensing: 

Propeller  thrust  can  also  be  considered  as  a  linear 
function  of  V0,  CO  ,  and  p  Just  as  Qp. 

T  =  K4*V0  +  K5*p  +  K6*  (96) 


where 


K4*=  *X 

^vo|  at  the  operating  point  (97) 


(98) 


k  *  =  a  I 

3     "^p  I  at  the  operating  point 


K6*  =  \Z 

^  I  at  the  operating  point  (99) 

*    *       -it- 
Values  for  K4  ,  K5  ,  and  Kg  are  found  exactly  as  were 


values  for  previous  stored  variables.   Again  using  Troost's 
curves: 

K  *  =  LT      KT(>n2D4  =  (..04)(1.024)ffi)2(l8.25)4 

4  AVq  "   (.2)nD      (.2) (10) (18,25) 

K4*  =  -1.325  x  103  lb-sec  ( 

ft 

K  *  _  AT  m       KT9n2p4  =  (08) (1.024) (3^)2(l8.25)4 

5  "  **  "   (.2)jD      (.2) (1.15) (18.25) 

KR*  a  6.00  x  103  lb 

5  ft   ' 

and 

*    ^T  _   K^n^D4  =  (.04)(1.024)(^)2(18.25)4 

6  *"    (.2)(n)(2TT)   (.2)dg)2T 


K**  =  2.47  x  104  lb-sec 
°  rad 

Incrementally,  equation  (96)  becomes 

T  =  K4*V£  +  K5V  +  H  '  <100) 

The  block  diagram  for  the  system  sensing  thrust  is  shown 
in  Figure  E-8. 

In  order  to  determine  the  effect  on  ^     of  thrust 
sensing  regulation  it  is  necessary  to  add  the  paths  that 
generated  cu     from  Vq,  p,  and  Q1.   This  is  shown  in  Figure 
E-9. 
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FIGURE  B-8 
System  With  Thrust  Sensing 
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FIGURE  E-9 

Block  Diagram  of  Complete  Thrust  Sensing  Regulation 

Showing  Generation  of 
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The  transfer  function  between  Vq(s)  and  u>  (s)  is 
more  complicated  in  this  case  due  to  the  presence  of  two 
paths  and  two  loops. 


*ilU) 


=    ATTrfl 


K4*GiK2 


L  »r  +  F^i7 


-] 


▼•(■)  1+K£GT  .  K6«k9«GT 

N(rps+i)      ^(TpS+iMYs+i) 


(101) 


Rearranging : 


u>' 


V0(b) 


-ki% 


la)  =  JE2 


-/El!        +  K2«K4«GT  \ 


^  +  ;2Tp     +  KS*Q- 


TpT 


WTt 


S+i 


K5»GT 

rv\TTp 


(102) 


-  K6*K2*GT 

M£rrrP 


Substituting  Values:  ( ~(,L\o  2..^6Xio*£l£25*io3^T 

6£|0     S 


I 


)^»*DCW)     ^6?  K-5^ 
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^  (s)  =     .01895  S  -    (70339  -  3190  GT) (103) 

Vo(sJ  S2+(2.09+1.7x105gt)S  +   (0. 545-7. 6xl03GT) 

Again  setting   J  =  0.4 

2  ^60n  =  2.09  +  1.7  x  105  GT 

£On  =  2.61  +  2.12  x  105  GT 

Solving  for  GT: 

(2.61  +  2.12  x  105GT)2  «  (.545  -  7.6  x  lO^Gj) 

4.5  x  1010GT2  +  4.31  x  105qt  +  5.26  =  0         (104) 


Equation  (104)  has  only  imaginary  roots.  In  fact 
setting  GT=0  in  equation  (103)  results  in 

^n  =  °«736  sec"1 
and  I  =  1.42       , 
as  well  as  a  gain  at  zero  frequency 


<o'(0) 


^>(0)   =  -  0.0621    ;   10  login  K 
V*(0)  V(0) 


=  -12.06 


which  checks  with  the  fixed  pitch  case  as  it  must. 

For  all  values  of  Gq,  greater  than  zero,  |  Increases 
and  zero  frequency  gain  also  Increases. 

For  negative  values  of  GT,  the  pitch  control  mechan- 
ism would  sense  a  decrease  in  thrust  and  respond  by  decreasing 
pitch  which  would  further  decrease  thrust  and  the  system 
would  settle  at  zero  pitch,  zero  thrust,  and  maximum  CO  . 

Assuming  a  value  of  Gm  =  +  10  J   £t  ,  the  unfortunate 

lb 

dilatorius  effect  of  the  thrust  sensing  can  be  shown. 


Equation  (103)  becomes 

W'(s)  =  +.01895.    .   _ 
Vo(s)    S2  +  3.79  S  +(.469) 


W\(b)   =  +.018955  ♦  .0658  (1Q5) 


Equation  (105)  is  plotted  in  Figure  E-6  for 
comparison  purposes. 
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